ARLINGTON  HALL  STATION 


ARLINGTON  12  VIRGINIA 


NOTICE:  WHEN  GOVERNMENT  OR  OTHER  DRAWINGS,  SPECIFICATIONS  CHI  OTHER  DATA 
'ffl’WD  FOR  ANY  PURPOSE  OTHER  THAN  IN  CONNECTION  WITH  A  DEFINITELY  RELATED 
GOVERNMENT  PROCUREMENT  OPERATION,  THE  U.  S.  GOVERNMENT  THEREBY  INCURS 
NO  RESPONSIBILITY,  NOR  ANY  OBLIGATION  WHATSOEVER;  AND  THE  FACT  THAT  THE 
GOVERNMENT  MAY  HAVE  FORMULATED,  FURNISHED,  OR  IN  ANY  WAY  SUPPLIED  THE 
SAID  DRAWINGS,  SPECIFICATIONS,  OR  OTHER  DATA  IS  NOT  TO  BE  REGARDED  BY 
IMPLICATION  OR  OTHERWISE  AS  IN  ANY  MANNER  LICENSING  THE  HOLDER  OR  ANY  OTHER 
PERSON  OR  CORPORATION,  OR  CONVEYING  ANY  RIGHTS  OR  PERMISSION  TO  MANUFACTURE, 
USE  OR  SELL  ANY  PATENTED  INVENTION  THAT  MAY  IN  ANY  WAY  BE  RELATED  THERETO. 


F 


I 


1 

I 


1  .  t 
\  '# 


>~ 

D_ 

CD 

CJ> 

loJ 


o 


Q  GO 


UNIVERSITY  OF  MINNESOTA 


SCHOOL  OF 

BUSINESS  ADMINISTRATION 


nu  co#»y 


kmtwrn  ♦* 

ASTI  A 

AtiiNOTON  hail  station 

AtUwCTOH  12  VIXOiHtA 
* 

tliSt 


A 

/• 


I  1 


JU'I  5W9 


June  1,  1959 


SOME  PROPERTIES  OP  THE  SPEARMAN  ESTIMATOR 

IN  BIOASSAY 

Byron  William  Brown,  Jr.  ^ 
Technical  Report  No.  6 


University  of  Minnesota 
Minneapolis,  Minnesota 

Reproduction  in  whole  or  in  part  is  permitted  for  any 
purpose  of  the  United  States  Government . 

^  Work  done  in  part  under  contract  Nonr2582(00) , 

Task  NR  042-200  of  the  Office  of  Naval  Research. 


iii 


TABLE  C!F  CONTENTS 


Page 


1.  INTRODUCTION  1 

1.1  An  Example  of  a  Quantal  Assay  1 

1.2  Terminology  2 

1.3  The  Experiment  and  the  Model  2 


2.  THE  SPEARMAN  ESTIMATOR  h 

2.1  Definition  of  the  Spearman  Estimator  1* 

2.2  The  Mean  and  Variance  of  the  Spearman 

Estimator  for  the  Finite  Experiment  6 

2.3  Comparison  of  the  Spearman  Estimator 

with  Parametric  Competitors  7 

2,U  Examples  of  Tolerance  Distributions  9 


3,  THE  INFINITE  EXPERT lENT:  THE  MEAN  AND  VARIANCE 

OF  Th£  SPEARMAN  ESTIMATOR  9 

3.1  The  Infinite  Experiment  9 

3.2  Mean  and  Variance  for  the  Spearman  Estimator 

for  the  Infinite  Experiment  II 


h.  EFFECT  OF  DOSE  MESH  LOCATION  ON  BIAS  AND  VARIANCE 


OF  THE  SPEARMAN  ESTIMATOR  IN  THE  INFINITE 
EXPERIMENT  13 

Uol  General  Discussion  13 

U,2  Bounds  on  the  Bias  of  the  Spearman 

Estimator  13 

U.2,1  Expression  for  the  Bias  111 

U.2.2  Bound  on  the  Bias  ill 

U.2,3  Tolerance  Distributions  for  Which 

the  Bound  is  Attained  Hi 

il,2,li  Bound  on  the  Bias  for  Unimodal 

Distributions  15> 

11.2.5  Unimodal  Tolerance  Distributions 

for  Which  the  Bound  is  Attained  17 

11.2.6  Bounds  on  the  Bias  for  Distributions 

with  Hi^er  Order  Derivatives  18 

11.2.7  Expressions  for  the  Bound  Involving 

Derivatives  of  the  Tolerance 
Distribution  21 

11.2.8  Computation  of  the  Bounds  on  the 

Bias  for  some  Tolerance  Distributions  23 
U,3  Effect  of  Dose  Mesh  Location  on  the 

Variance  of  the  Spearman  Estimator  2ii 

U.3*l  Bound  on  |V(x|xo)  -  V|  2lt 

ii.3,2  A  Tolerance  Distribution  for  Which 

I  V  (x  }  x  )  -Vi  Approximates  the 
Bound  Arbitrarily  Closely  2? 


iv 


Page 

4.3.3  Bounds  on  j  V(x~|x0)  -  V  I  in 
Tertce  of  Derivatives 

28 

4.3*4  Computation  of  Bounds  on 
|V(x  |xQ)  -  V  )  for  some 

Tolerance  Distributions  29 

5«  RANDOM  LOCATION  OF  THE  DOSE  MESH  IN  THE 

INFINITE  EXPERIMENT  30 

5*1  Introduction  30 

$•2  Unbiasedness  of  the  Spearman  Estimator  30 

5*3  The  Mean  Square  Error  of  the  Spearman 

Estimator  32 

5*4  Values  of  V  for  Several  Tolerance 

Distributions  34 

6,  LARGE  SAMPLE  PROPERTIES  OF  THE  SPEARMAN 

ESTIMATOR  FOR  THE  INFINITE  EXPERIMENT  '  34 

6.1  Large  Sample  Experiments  34 

6.2  Optimum  Choice  of  n  and  d  for  the 

Spearman  Estimator  36 

6.3  Large  Sample  Properties  of  the  Spearman 

Estimator  41 

7.  URGE  SAMPLE  EFFICIENCY  OF  THE  SPEARMAN 

ESTIMATOR  FOR  THE  INFINITE  EXPERIMENT  42 

7.1  Previous  Comparisons  of  the  Spearman 

Estimator  with  the  Maximum  Likelihood 
Estimator  42 

7.2  Definition  of  Asymptotic  Efficiency  43 

7.3  The  Spearman  Efficiency  for  Several 

Tolerance  Distributions  44 

7.3*1  Logistic  44 

7.3*2  Normal  44 

7.3.3  Angular  44 

7.3.4  Uniform  45 

7.3.5  One  Particle  45 

7.3.6  Algebraic  45 

7.4  The  Logistic  Tolerance  Distribution  and 

the  Spearman  Estimator  47 

7.5  Distributions  with  Efficiency  of  the 

Spearman  Estimator  Close  to  Zero  $0 

7.6  Two  Parameter  Families  of  Tolerance 

Distributions  52 

8.  THE  SPEARMAN  ESTIMATOR  FOR  FINITE  EXPERIMENTS  54 

8.1  Finite  Experiments  54 

8. 2  The  Spearman  Estimator  54 

8.3  Information  for  the  Finite  Experiment  55 

8.4  Efficiency  for  the  Spearman  Estimator  for  the 

Finite  Experiment,  Scale  Parameter  Known  60 

8.5  Efficiency  for  the  Spearman  Estimator  for  the 

Finite  Experiment,  Scale  Parameter  Unknown  60 


V 


Page 

9.  REGULAR  EEST  ASYMPTOTIC  ALII  NORMAL  ESTIMATES 

WITH  THE  WRONG  MODEL  62 

9*1  General  Discussion  6*2 

9*2  One  Level  Experiment  67 

9.3  Two  Level  Experiment  73 

9.U  The  2k+l  Level  Experiment  77 

9.$  Suianary  of  Wrong  Model  Investigation  82 

BIBLIOGRAPHT  83 

APPENDIX  Is  THE  DISTRIBUTION  FUNCTION  AND 

INFORMATION  FOR  THE  INFINITE  EXPERIMENT  86 

1*1  Distribution  Function  for  the  Infinite 

Experiment  86 

I #2  .  Information  for  the  Infinite  Experiment, 

Scale  Parameter  Known  86 

I. 3  Information  for  the  Infinite  Experiment, 

Scale  Parameter  Unknown  87 

APPENDIX  II s  BERNOULLI  PERIODIC  FUNCTIONS  AND  THE 

EULER-MAOLAURIN  FORMULAE  90 

II.  1  The  Bernoulli  Periodic  Functions  90 

II.  2  The  Euler-MacLaurin  Formulae  91 

APPENDIX  III:  ESTIMATION  OF  TiE  VARIANCE  OF  THE 

SPEARMAN  ESTIMATOR  __  93 

III. l  An  Estimator  for  V(3c)  93 

III ,2  An  Alternative  Estimator  of  V  Based  on 

the  Second  Moment  of  F  9I4 


i 


Vi 


GLOSSARY 

1.  The  Tolerance  Distribution 


x 

F(x) 
f(x) 
f (n)(x) 


a  dose  level 

•a  tolerance  distribution 

F»(x) 

the  nth  derivative  of  f 

the  mean  of  the  tolerance  distribution 


m 


L.5o 


a 

R 


m 


2.  The  Sxperimem 


the  mode  of  F 

the  median  of  F 

the  standard  deviation  of  F 

the  distance  between  the  20th  and  80th 
percentiles  of  F 

f(xm),  the  maxinrum  of  f(x) . 


x. 

i 


n. 

x 


n 


the  ith  dose  level 


the  number  of  subjects  tested  at  x^ 

the  number  of  subjects  tested  at  each  x^ 
when  •'  he  sample  sizes  are  all  equal 

the  common  distance  between  dose  levels: 
d=x.  -  ,-x. 

X+l  X 


d/<7 

d/R 


x_ 


k 


the  "middle "  dose  level  (x.*x  +id, 

x  o 

the  number  of  dose  levels  on  each  side  of 
xq  for  the  finite  experiment 

the  distance  on  each  side  of  x  covered 

o 

by  dose  levels  in  the  finite  experiment: 
a*kd 


3. 


L 


a 


▼u 


N  :  the  total  number  of  subjects  tested  in 

the  finite  experiment :  M«n(2k+1) 


3.  Observations  and  Estimators 


ri 

• 

• 

the  observed  number  of 
at  the  ith  dose  level. 

subjects  responding 
xi 

pi 

• 

• 

the  propertidn  of 

at  the  ith  d03e  level: 

subjects  responding 

wS. 

X 

• 

the  Spearman  estimator 

ii.  Information  and  Characteristics  of  Estimators 


E  (  ) 

V  (  ) 

MSE  (  ) 

B  (  ) 


denotes  expectation  when  a  random  variable 
appears  in  the  brackets 

denotes  variance 

denotes  mean  square  error 

denotes  bias 


subscript 
a  or  k 


E  (  ) 

x 

o 


subscript  A 


7 


denotes  finite  experiment  (doses  x  «x  +ld, 
i*Q,  ±1,  ±2,....±k).  If  no  subscrlpt°a  or 
k  appears,  the  experiment  is  infinite 

(doses  x.»x  +id,  i-0,  ±1,  ±2, . ).  If 

no  conditional  notation  concerning  x 
appears,  x  is  taken  to  be  randomly  Bhosen 
on  the  interval  (o,d). 

expectation  with  respect  to  x  over  the 
interval  (0,d)  0 

denotes  an  asymptotic  moment  (see  (6.3) 
and  (9.1)) 

the  value  j- j  F(l-F)dx*  7  is  shown  in 

section  ...3  to  approximate  V(x)  and  is 
defined  in  section  6,3  to  be  the  asymptotic 
variance  of  x  for  the  infinite  experiment 

information  for  scale  parameter  known 

element  in  the  inverse  information  matrix 
corresponding  to  p.  when  the  scale  parameter 
is  unknown 
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E  :  denotes  efficiency  relative  to  I  if  no 

estimrtor  follows  in  brackets.  E(x  ) 
denotes  efficiency  conditional  on  x°. 

11  11 
E  :  denotes  efficiency  relative  to  I 


1.  INTRODUCTION 

1*1  An  Example  of  a  Quantal  Assay: 

In  certain  experimental  situations  it  is  impossible  to 
measure  the  variable  on  each  experimental  unit  directly,  but 
it  is  possible  to  fix'*'  a  number  for  each  unit  and  then  deter¬ 
mine  simply  whether  or  not  the  experimental  unit  has  a  measure¬ 
ment  greater  than  this  number.  This  type  of  experimental  situ¬ 
ation  is  found  in  a  variety  of  fields  of  biological  investiga¬ 
tion.  The  following  example  from  hormone  assay  illustrates 
the  nature  of  the  problem(lS). 

Some  estrogenic  preparations  are  extracts  from  the  urine 
of  pregnant  mares.  These  preparations  are  mixtures  of  several 
estrogens.  The  estrogenic  strength  of  such  extracts  can  not 
be  measured  analytically  in  a  satisfactory  way.  It  is  known, 
however,  that  if  sufficient  estrogenic  substance  is  given  to 
immature  or  spayed  female  mice,  they  will  show  cornification 
of  the  vagina.  If  a  fixed  dose  of  the  preparation  is  given 
to  a  test  animal  and  cornification  is  observed,  all  that  is 
known  is  that  the  dose  administered  was  at  least  as  great  as 
the  ’tolerance*  of  the  mouse  to  this  preparation.  If  no 
cornification  is  observed,  there  is  no  way  of  knowing  how 

^In  some  applications  this  number  can  be  fixed  only  with 
appreciable  error.  This  case  is  not  considered  in  this  paper, 
but  has  been  discussed  by  Haley  (21)  for  certain  parametric 
formulations . 
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such  higher  the  dose  would  have  had  to  be  to  induce  the 
response.  It  is  desirable  to  estimate  the  mean  tolerance 
of  a  group  of  test  animals  to  the  preparation  on  the  basis 
of  such  data.  The  mean  tolerance  and  the  strength  are  in¬ 
versely  related.  Other  examples  of  this  assay  situation  can 
be  found  in  insecticide  research  (7),  vitamin  research  (20), 
vaccine  screening  for  safety  (25),  and  toxicity  evaluation 
of  various  chemicals  (27).  The  quantal  assay  situation  can 
also  be  found  in  industry  in  munitions  testing  (14)  and 
reliability  testing  (30)  among  other  applications. 

1.2  Terminology 

The  variable  under  investigation  will  be  called  the  dose. 
The  dose  may  be  a  direct  measure  of  the  stimulus  (e„g.  con¬ 
centration  of  an  injection)  or  it  may  be  some  transformation 
of  this  direct  measurement  (comaonly  the  log  of  the  measure¬ 
ment),  The  experimental  units  will  be  tested  at  various 
doses.  The  observation  on  each  unit  will  be  either  a  response 
or  a  no-response,  Jlll-or-none  responses  are  called  quanta  1 
responses  in  biological  experiments.  The  experiment  is  av 
quantal  assay.  The  probability  of  response  depends  on  the 
dose.  The  function  relating  the  probability  of  response  to 
the  dose  level  Is  the  tolerance  distribution  or  dose  response 
function. 

1.5  The  Experiment  and  Model 

The  usual  quantal  assay  is  done  in  the  following  manner. 
k  set  of  dose  levels  is  chosen.  The  test  subjects  to  be  uaed 
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in  the  experiment  are  randomly  allocated  to  the  dose  levels. 

The  number  of  responses  among  the  subjects  in  each  dose  level 
group  is  used  as  the  basis  for  inference. 

The  test  animals  should  be  randomly  selected  from  a  well- 
defined  population  of  subjects.  The  dose  response  function 
is  descriptive  of  this  population.  The  stability  of  the  dose 
response  function  over  time  must  be  investigated  (3). 

The  notation  used  in  this  paper  will  be  as  follows: 

x  denotes  a  dose  level 

F(x)  denotes  the  dose  response  function,  i.e.  the 
expected  proportion  of  responses  at  dose  x 
denotes  the  i^1  dose  level  used  in  the  experiment, 

1*0,  ±1,  ±2,... 

denotes  the  number  of  subjects  tested  with  dose  x^ 
denotes  the  number  of  subjects  responding  among  the 
subjects  receiving  dose 

denotes  the  proportion  of  subjects  responding  among 
the  group  of  subjects  receiving  dose  x.,  p.  *  i/ 

jL  X 

The  assumptions  concerning  the  dose  response  function 
and  the  observations  are: 

F(x)  is  a  distribution  function 

F(x)  has  a  first  moment,  p,  called  the  mean  tolerance 

The  observations  on  the  subjects  are  mutually  independent 
The  observations  are  the  dichotomous  quantal  response  vari¬ 
ables.  A  set  of  sufficient  statistics  for  the  experiment  oawdsta 


If 


\ 


{ 


u 


of  the  numbers  of  responses,  r^,  among  the  subjects  tested  at 
each  dose  level.  The  r^  are  mutually  independently  distributed 
binomial  variables  with  means  niF(xi),  i»0,  ±1,  ±,..0 

The  primary  problem  in  this  paper  is  the  estimation  of  p 
on  the  basis  of  the  experiment  described  above.  Some  bio- 
statistical  writers  (U,  5,  6,  8,  10,  15,  16,  17,  18)  recommend 
parametric  estimators  for  this  problem,  i.e,  estimators  which 
necessitate  the  specification  of  a  functional  form  for  the 
dose  response  function.  One  function  used  frequently  is  the 
normal  tolerance  distribution: 

,  ££ 

F(x;n,o2)  -  -L-  \  a  -  t2/„  dt  (1.1) 

.IT?  e  2 

J 

This  paper  presents  an  evaluation  of  the  Spearman  estimator, 
a  nonparametric  estimator  of  p. 

2.  THE  SPEARMAN  ESTIMATOR 

2.1  Definition  of  the  Spearman  Estimator 

The  estimator  to  be  discussed  in  this  paper  was  des¬ 
cribed  by  Spearman  (28)  in  1908.  He  gives  credit  for  the 
idea  to  the  German  psycho-physiologist,  Muller.  The  esti¬ 
mator  was  described  again  by  Karber  (23)  in  1931  and  is 
occasionally  referred  to  as  the  Spearman-K&rber  estimator. 

Spearman  defined  the  e  stimator  for  regularly  spaced 
dose  levels,  x^*XQ+id,  i»0,  ±1,  ±2,  and  equal  num¬ 

bers  of  subjects  tested  at  a  finite  number  of  dose  levels, 
say  n^*n  for  i*0,  ±1,  ±2,  ...±k.  The  estimator  is 
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k-1 

2 

-k 


(x,  -  72) 


(pi+l 


-Pi). 


x  is  analogous  to  a  grouped  mean  for  continuous  data. 

Spearman  computed  this  estimator  only  when  p_^  *  0 
and  p^  *  1,  In  practice  the  following  modification  is  used: 


k-1 


p-k(*-k-  *  5k  <v  /2)  <pi+rPi) +  /2~>  (2.i) 


Thus  any  estimate  of  probability,  p_^,  below  the  lowest  do3e 
level  (x_^)  for  which  n^/  0  is  assigned  to  the  point  half  a 
dose  interval  below  the  lowest  levelj  and  the  estimate  of 
probability,  1-p^,  above  is  handled  similarly. 

Armitage  and  Allen  (2)  extended  Spearman's  definition 
to  unequally  spaced  dose  levels,  x^: 


k-1  x.+x... 

5  .  Z  (  (Pitl-Pi) 

-k 


(2.2) 


This  definition  can  be  modified  to  allow  for  estimates  of 
probability  below  and  above  x^. 

Irwin  (22)  and  Finney  (15,  16)  discussed  an  experiment 
in  which  subjects  would  be  tested  at  an  infinity  of  dose 
levels.  This  experiment  called  for  n^*n  and  x^*xQ+id, 
i»0,  ±1,  ±2,...,  with  xq  and  d  chosen  arbitrarily.  (See 
Appendix  I  for  a  discussion  of  the  resulting  infinite 
sample  space,) 

If  n./o  for  all  i  and  x.  are  chosen  so  that  x.-*«  as 
x  i  i 

i — and  x. — >  -*•  as  i — }  ,  then  the  estimator  can  be 

x 

defined  as  the  limit  of  x  defined  in  (2.2)  as  k-^-*  , 


x  »  lira 


6 


k-1 

Z 

-k 


(2.3) 


M  ( - 


The  following  experimental  designs  will  be  considered: 

(a)  •  xq  +  id  with  n^  -  n  for  i  *  0,  ±1,  ±2,  ,  ,.±k 

and  ni  *  G  for  i  «  ±(k+l),  ±(k+2),.,..  (b)  xi  ■  xq  +  id 
with  n^  *  n  for  i  *  0,  ±1,  ±2........  The  first  experi¬ 
ment  will  be  referred  to  as  the  finite  experiment  and  the 
second  experiment  will  be  referred  to  as  the  infinite  exper- 
iment. 


Notes  The  Spearman  estimator  for  the  finite  experi¬ 
ment  (2.1)  can  be  expressed  in  several  ways: 

5  \  *  %  -  d  Pi  <2 M 

,  k  0 

*  *  Xo  +  %  +  d  2  qi  ‘  d  2  pi  (2.5) 

1  -k 

qi  *  1  -  pi 

2.2  The  Mean  and  Variance  of  the  Spearman  Estimator  for  the  Finite 
Experiment 

The  exact  mean  and  variance  of  the  Spearman  estimator 
for  the  finite  experiment  and  fixed  ?:  are: 


Vx|x„) 


<*-k-d/2  >p.k+kM  vd/2  >  i-Fi>+<  vd4  >  (1-Fk> 


where  F.  =  F(x. ) 

l  i 


(2.6) 


k 

2  F.(1-FJ 
-k  1  1 


(2.7) 


As  an  illustration  (Table  2.1),  the  bias  and  the  variance  of 
the  Spearman  estimator  have  been  computed  for  a  normal  tol¬ 
erance  distributi’v  '1.1).  The  experimental  design  consists 
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of  five  dose  levels,  two  standard  deviation  units  apart, 
with  n  subjects  at  each  of  the  five  levels,  '"‘he  bias  and 
the  variance  of  the  estimator  depend  on  th'"  .ocation  of 
the  dose  mesh  relative  to  the  mean  of  the  olerance  dis¬ 
tribution,  Therefore  the  bias  and  variance  were  computed 
for  several  locations  of  the  dose  mesh.  The  location  of 
the  dose  mesh  is  indicated  by  the  distance  of  the  middle 
dose,  x  .  from  the  mean  of  the  tolerance  distribution,  u, 

'  O' 

in  standard  deviation  units. 

When  the  mean  is  within  the  interval  spanned  by  the  dose 
levels,  the  fluctuations  in  the  bias  and  the  variance  as 
functions  of  the  location  of  the  dose  mesh  are  negligible. 
When  the  dose  mesh  fails  to  cover  the  mean  the  bias  becomes 
large  and  the  variance  goes  to  zero.  The  mean  square  errors 
for  the  case  of  n  equal  to  ten  and  to  one  hundred  are  also 
shown  in  Table  2.1. 

2.3  Comparison  of  the  Spearman  Estimator  with  Parametric 
Competitors 

It  is  clear  from  the  above  introduction  to  the  Spearman 
estimator  that  it  has  certain  advantages  over  its  parametric 
competitors : 

a)  The  Spearman  estimator  is  simple  in  concept,  being  just 
the  mean  of  a  histogram  reconstructed  from  the  quanta! 
data, 

b)  The  Spearman  estimator  is  simple  to  conpute,  It  involves 
only  the  sum  of  the  observed  proportion*  (2.1;),  The 
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Tablo  2.1 


The  Bias,,  Variance  and  Mf  an  Square  Error  of  the  Spearman  Estimator 
(2.1)  of  the  Mean  of  a  Dermal  Tolerance  Distribution,  Using  Five 


Dose 

Levels  Spaced  2a  Apart 

X  -a 
o  r 

Bias 

nV*  i  V 

. -y—  "■ 

nMSE  (x 

j  X  )/o2 

e 

a 

2 

O 

n*10 

n-100 

0 

0 

1.178 

1.178 

1.178 

.2 

-.003 

1.166 

1.168 

1.169 

•  ii 

-.00U 

1.1UU 

l.liiU 

I.U16 

.6 

-.005 

1.11U 

l.llii 

1.117 

.8 

-.003 

1.089 

1.089 

1.090 

1.0 

0 

1.078 

1.078 

1.078 

1.2 

.003 

1.089 

1.089 

1.090 

l.U 

.005 

l.llii 

1.111* 

1.117 

1.6 

.005 

l.liiU 

l.liUi 

1.1U7 

'  1.8 

.003 

1.168 

1.168 

1.169 

2.0 

0 

1.178 

1.178 

1.178 

2.2 

-.00  3 

1.168 

1.168 

1.169 

2.U 

-.OOii 

1.1U3 

1.1U3 

1.1U5 

2.6 

-.OOLi 

1.112 

1.112 

l.llii 

2.8 

-.001 

1.086 

1.086 

1.086 

3.0 

.003 

1.073 

1.073 

1.07U 

3.2 

.008 

1.078 

1.079 

1.081* 

3.U 

.oil* 

1.095 

1.097 

1.115 

3.6 

.021 

1.112 

1.116 

1.156 

3.8 

.031 

l.llii 

1.12U 

1.210 

h.o 

.oU6 

1.089 

1.110 

1.301 

h.2 

.069 

1.030 

1.078 

1.506 

h.h 

.106 

.936 

1.0U8 

2.060 

U.6 

.158 

.815 

1.065 

3.311 

U,8 

,229 

.678 

1.202 

5.922 

5.0 

.320 

.539 

1.563 

10.779 

5-2 

.U32 

.Uio 

2.276 

5.U 

.563 

.299 

3.U69 

5.6 

.710 

.208 

5.2U9 

5.8 

.872 

.139 

7.7U3 

6.0 

I.0I46 

.089 

II.030 

6.2 

1.228 

.055 

6.U 

1.U16 

.032 

6.6 

1.609 

.019 
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parametric  estimation  procedures  ordinarily  involve 
either  an  iterative  solution  or  a  weighted  regreael on 
solution. 

c)  The  exact  mean  and  variance  of  the  Spearman  estimator 
ars  easily  obtained  for  any  size  sample.  Therefore 
experimental  design  investigations  are  readily  done. 

In  contrast,  only  asymptotic  theory  is  available  for 
the  parametric  estimators. 

d)  The  Spearman  estimator  is  nonparametric  in  that  no 
functional  form  need  be  assigned  F(x)  in  order  to  com¬ 
pute  the  value  of  the  estimate  from  the  data.  How- 
important  this  point  is  depends  on  the  robustness  of 
the  parametric  estimators.  Some  relevant  results  are 
presented  later  in  this  paper  (section  9). 

2.ii  Examples  of  Tolerance  Distributions 

The  tolerance  distributions  in  Table  2.2  serve  as  illu¬ 
strations  throughout  this  paper.  The  first  four  distribu¬ 
tions  are  used  in  practice.  The  remaining 

distributions  are  to  illustrate  specific  points  (see 
sections  U.2,  U.3,  5.U,  7.3,  and  7.6). 

3.  THE  INFINITE  EXPEFJ32NT :  THE  MEAN  AND  VARIANCE  OF  THE  SPEAR¬ 
MAN  ESTIMATOR 

3.1  The  Infinite  Experiment 

Irwin’s  (22)  and  Finney’s  (l£,  16)  concept  of  an  iiifinite 
experiment  (section  2.1)  makes  possible  a  mathematical  dis¬ 
cussion  of  the  effect  of  the  location  of  the  dose  mesh  on 
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Table  2.2 

Tolerance  Distribution* 


NAME 

FuncnonAL  form 

VARIANCE 

1.  Normal 

x  p  P(x-tl)  e-t2/2 

7*2 

dt;  -®<x<od 

1 

I* 

2.  Logistic 

-00 <x<® 

«* 

— 5 

302 

3.  Angular 

•ln2[p(x-(i)+*/4]j 

-*A^(*-y)<*A 

x2-8 

Kb  Cbe-Partid* 

X>0 

16? 

y2 

5 •  Uniform 

3(x-n)+  ^2  ; 

-1//2<0(x-m)  <  'l/2 

l 

180s 

6,  Algebraic 

,  -8 

1-x  ; 

X>1,  SOI 

s 

(s-2)(s-l) 

7.  Student* s 

r  (*-y) 

k.  J 

-oo  dt 

-®<x<CD,e>0 

2e+l 

uttr 

/-  t2  vl4e 
(i 

o 

Hote  1:  For  the  algebraic  distribution  the  mean  is  s-I 

• 

For 

the  other  distributions 

the  mean  is  ix. 

Note  2:  For 

the  first  five  distributions  0  must  be  positive  and 

the 

variance  exists  for  all  0.  For  the  last  two 

distributions  the  variance  exists  only  if  s>2  and  e>| 
respectively. 


u 


the  bias  and  the  variance  of  the  estimator  without  discussing 
the  possibilities  of  grossly  misplacing  the  whole  set  of  dose 
levels  relative  to  p  •  The  concept  of  the  infinite  experi¬ 
ment  also  facilitates  the  development  of  large  sample  defini¬ 
tions  of  mean  square  error  and  efficiency  in  later  sections 
(6  and  7).  The  investigation  of  the  infinite  experiment 
has  practical  importance  since  it  is  shown  in  Appendix  I 
that  the  information  for  the  infinite  experiment  is  essential¬ 
ly  the  same  as  that  for  the  corresponding  finite  experiment 
covering  ''most"  of  the  range  of  F(x). 

3.2  Mean  and  Variance  of  the  Spearman  Estimator  for  the  Infinite 
Experiment 

Lemmas  3.2,1,  3.2.2,  and  3.2,3  establish  conditions 
under  which  the  Spearman  estimator  has  the  following  mean 
and  variance: 


E(x|xQ)  -  2(xi+d/2)(Fi+1-Fi) 

— oc 

(3.D 

00 

v(xixo}  "  r  2Fi<1"pi) 

(3.2) 

Using  (2,5),  the  estimator  is: 

x  “  lin  {  x0+d/o+d^“d^  Pi] 
k  — <•  c  1  1  -k  1 ' 


(3.3) 


Similarly 


mOb  «■*•© 


.1*) 


Lemma  3.2.1:  If  F  has  a  first  moment,  | x,  then  the  series 

X0+/2+d2(1"Fi)-d2Fi  (3.S) 


converges  to  a  finite  value. 
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Proof:  The  lemma  is  established  if  dZ(l-F^)  and  dZF^ 

are  shown  to  converge  to  finite  values.  Consider  the 

remainder  for  dZ(l-F. )  : 

1  1 


d  I  (1- 
k+1 


-V  *  j 


(l-F)dx 


Interchanging  the  order  of  integration  on  the  right: 

m  #  m 

d  2  (1-F. )  -  (  tdF(t) 
k+1  1  ) 

Tc 

If  F  has  a  first  moment,  the  integral  on  the  right  goes  to 

zero  as  x.  becomes  infinite.  Therefore  dZ(l-F. )  is  finite. 

K  o  1  x 

Similarly  dZF^  can  be  shown  to  be  finite.  Q.E.D. 

Lemma  3.2.2:  If  F  has  a  first  moment,  then  the  series 


-  ZF, (1-F.) 
n  iv 


(3.6) 


converges. 

Proof:  In  the  proof  of  Lemma  ^.2.1  it  was  established  that 

*  ^ 

Z(l-F^)  is  finite.  Since  0  ■  F^l-F.^)  -  1-F^  the  series 

*>  o 

ZF.(1-F. )  also  is  finite.  Similarly  ZF, (1-F.)  is  finite, 

i  X  * 0*  ^  X 

d2  * 

and,  therefore,  —  ZF^(l-F^)  is  finite.  Q.E.D. 

lemma  3.2.3:  If  F  has  a  first  moment,  the  Spearman  estimator  (3.3) 
for  the  infinite  experiment  converges  with  probability  one 
to  a  random  variable  with  mean  and  variance  given  by  (3.1) 
and  (3.2). 

Proof;  Leninas  3.2.1  and  3*2.2  prove  the  convergence  of  the 
series  in  (3.5)  and  (3.6).  Theorem  2,3  in  chapter  III  of 
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Doob  (12)  establishes  the  sufficiency  of  the  convergence 
of  these  two  series  for  the  convergence  of  x  with  probabil¬ 
ity  one,  and  shows  that  the  expected  value  and  variance  of 
x  are  given  respectively  by  (3.5)  and  (3.6),  or,  using  (3.U), 
by  (3.1)  and  (3.2).  Q.E.D. 

The  variance  of  the  Spearman  estimator  is  not  zero  even 
though  the  number  of  observations  is  infinite.  Similarly, 
it  is  shown  in  Appendix  I  that  the  information  contained  in 
the  infinite  experiment  is  finite  in  the  common  parametric 
formulations. 

EFFECT  OF  DOSE  MESH  LOCATION  ON  BIAS  AND  VARIANCE  OF  THE 
SPEARMAN  ESTIMATOR  IN  THE  INFINITE  EXPERIMENT 

2.  General  Discussion 

In  the  infinite  experiment  the  expected  value  and  vari¬ 
ance  of  the  estimator  depend,  in  general,  on  the  doses, 

x^*xQ+id,  1*0,  tl,  ±2, .  For  a  particular  spacing,  d,  the 

expected  value  and  variance  will  be  simply  functions  of  xq 
with  period  d. 

Finney  (15,16)  has  computed  the  bias  and  variance  for 
normal  and  logistic  tolerance  distributions.  Since  the  Spear¬ 
man  estimator  is  nonpar ametric,  it  is  desirable  to  have  infor¬ 
mation  on  the  bias  over  a  wide  class  of  distribution  functions. 
It  is  possible  to  find  bounds  for  the  bias,  distributions  that 
maximize  the  bias,  and  conditions  on  F(x)  that  limit  the  bias. 
Bounds  for  the  fluctuation  in  the  variance  can  also  be  obtained. 

„2  Bounds  on  the  Bias  of  the  Spearman  Estimator 


4,2*1  Expression  for  the  Bi as 

The  bias  of  an  estimator  is  the  difference  between  its 
expected  value  and  the  parameter  estimated.  Denote  the 
bias  for  the  Spearman  estimator  for  the  infinite 
experiment,  conditional  on  x  ,  by  B(x|xq).  Using  (3.1): 

B(x|xq)  -  E(xixo)-p 

.j(Xi+d/2)(riti-F.)-  f  xdF(x) 

««o 

oo  ^  xl+l 

*  2(xi+d/2)(F.+1-Fi)-2  f  xdF(x) 

_oo  — oo  J 

X. 

1 


2(xi+/rCi)(Fi+i-Fi) 


(4.1) 


i+1 


where 


c  .  =  — 

‘  s 


J  xdF(x) 

_ _  •  (4.2) 


Xi+1 


dP(x) 


4.2,2.  Bound  on  the  Bias 

From  (4.1)  it  follows  that: 
Lemma  4.2.2 


|B(x}xo)|^;.txi+d/2- 

»QD 


<  d 

*  ? 


(4.3) 


4.2.3,  Tolerance  Distributions  ,br  which  the  Pound  is  Attained 
The  bound  given  in  (4.3)  is  attained,  e.g.  the  one 
point  distribution.  In  this  case  when  one  of  the  dose  levels 
coincides  with  the  mass  point  of  the  distribution,  the  true 
mean  of  the  distribution  equals  the  dose  level,  but  the  esti¬ 
mator  has  an  expected  value  d/2  units  below  this  dose  level 
(at  the  midpoint  of  the  dose  interval  showing  the  probability 
increment).  This  example  indicates  that  for  any  distribution 
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function,  if  d  is  large  relative  to  the  dispersion  of  F(x), 
the  bias  can  be  approximately  d/2» 

Even  if  one  excludes  the  one  point  distribution  a  bias 
of  /g  is  attainable.  Consider  the  class  of  discrete  dis¬ 
tributions  with  mass  points  on  a  lattice  with  spacing  D. 

Then  if  the  dose  mesh  has  d«D/m,  m  a  positive  integer,  and 
if  the  dose  mesh  is  located  so  that  the  mass  points  coincide 
with  dose  levels,  then  each  (h.2)  has  the  value  x^  and 

B(x|xo)  -  f  ZCF^-Fp  -  | 

mmqo 

4.2.4.  Bound  on  the  Bias  for  Unimodal  Distributions 

The  situations  discussed  in  paragraph  4.2. 3  do  not  occur 
often  in  practice.  In  this  paragraph  F(x)  is  restricted  to 
functions  with  the  usual  properties  possessed  by  tolerance 
distributions • 

Lemma  4.2.4.  If  F(x)  has  a  unimodal  density  with  maximum  ordinate, 
f  ,  then 

nr  2 

lB(Sixo>!  *  r 

(Unimodality  means  the  densi.,  ,  f(x),  is  non-decreasing  for 
x  less  than  the  mode,  x  ,  ano.  non-increasing  for  x  greater 
than  the  mode.) 

Proof:  The  bound  can  be  obtained  by  examining  the  terms  in 
the  expression  for  the  bias  (U.l) , 

Let  I  «(x.,x.  .)  be  an  interval  located  above  the  modal 

i  i'  l+l 

value  x  ,  so  that  f  (x)  is  non-increasing  in  I, .  Then,  using 
m'  i 

c  as  defined  in  (4.2), 

X 
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<  ..  .d 


c  ^  *  x^+  /g  no  matter  what  the  density  f(x)  is  in  1^. 

For  the  given  probability  F^^-F^  associated  with  1^ 
the  minimum  value  possible  for  c  ^  is  attained  for  the 
density,  g: 

g(x)  =  f .  x.  *  x  »  x.+Rd 

&  i  x  x 

sM  -  fi+1  xi+P.d<  x<  xi+1 


R  is  determined  by: 


Xi+1 


g(x)dx  -  Fi+1-Fi  , 


i.e.  Rd(Vfi,l)*dfi+l  *  Vl-Fi  * 

For  the  minimizing  density,  g(x),  the  value  of  c.^  becomes: 

lg  F.  ,-F. 

x+1  l 

Since  for  the  density  f(x)  is  bounded  below  by  and 

above  by  x^+  /  ,  the  bias  term  for  1^  satisfies 


o  i  (x>d/2-  ei)(Fltl-?1)  S  (x^%-  =ig)(Fitl-F.) 


The  term  on  the  right  can  be  evaluated  in  terms  of  d,  R, 
and  the  values  of  f(x)  to  o':  tain: 

0  s  <Vd/2-  5  £  Rd-HXq-W 

Since  0  *  R  *  1,  this  inequality  can  be  relaxed  to  obtain: 
0  i  (xi+<i/2"  °i)(Fi+l-5'‘)  £  T-  (U-6> 


8  '  i  i+1' 


17 


Similarly,  for  intervals  belov  the  node,  the  following 
inequality  is  obtained: 

•  s  <*ui*1/2-ci)<WFi>‘  0 


Also,  for  the  interval  containing  the  mode,  say  (x^x^), 
the  inequality  obtained  is: 


sWfo> s  +d4-c0Hvpo> s  r  <Vfi>  <4-8> 


Then  (4*6) ,  (4.7),  and  (4.8)  combine  to  give  a  bound  on  the 
sum  of  the  contributions  to  the  bias  from  all  of  the  intervals: 


|B(x|x  )|  S  max  J 


CO 


(4.9) 


i-  (f  -f.)  ♦£  d‘  ,  , 

°  m  1  i  r  (fi*W 


r  *L  r  (fi+rfi) 

- CD 


oo  0 

Then,  since  £  (f,-f, and  £  (f.  .-f.)«f  ,  the  bound  is 

i  i+1  1  _  i+1  x'  o 

i  -00 


|b(;i*0)I  *  §-  fm 

Q.  E.  D. 

4.?. 5.  Uniaodal  Tolerance  Distributions  for  Which  the  Bound  is 
Attained 

The  bound  (4.4)  is  attained  for  certain  unimodal  densities. 
Consider  the  following  example:  Let  f(x)«l  for  0*xSl,  f(x)»0 
otherwise.  Let  the  dose  interval,  d,  be  given  by: 


d  m  - - -  (N  a  positive  integer) 

» *VB 

Let  xQ«0.  The  the  dose  levels  will  be  given  by: 
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N 

777 

»t+  f  2 


N+l 

^4 


*  •  • 


for  all  intervals  except  the  one  containing  the  point,  x*l, 
f(x)  is  uniform  and  c.^x^+'Vg  so  that  the  contributions  to 
the  bias  are  zero.  For  the  interval  containing  x-1,  i.e. 
the  interval 

/  N  N+l 

N+V2  *  N +V2 


the  contribution  to  the  bias  can  be  computed  as  follows: 
N 


1  N+l  N  N+1/ 

=  -  +  -  ( — -  -  — - )  -  —— 

N  N+/0  2  N+1/0  N+V9  N+A 


c  = 


h 


d  /  , 
XN+  /2=1 


F„  ,-F 


Hi  « 


2(:i+  /2 ) 


B(x|xJ» 


rN+1/l 

N+1/. 


-1 


1  11 

vN+V2)  (N+V2)2  8 


This  is  the  bound  given  in  (u.U)  since  d*  — 


N+1/. 


and  f  *1. 
m 


The  uniform  distribution  has  the  properties  of  symmetiy 
and  only  two  points  of  inflection.  Therefore,  these  pro¬ 
perties  do  not  lead  to  a  stronger  bound  on  the  bias  of  the 
Spearman  estimator, 

U.2.6.  Bounds  on  the  Bias  for  Distributions  io  Tern*  of 
Derivatives 

The  example  of  tne  uniform  distribution  in  paragraph  lx. 2.$ 
suggests  the  contributions  to  the  bias  of  the  Spearman  esti- 
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mator  come  from  discontinuities  and  rapid  rates  of  change 
in  the  density  function.  3ounds  on  the  bias  can  be  tight¬ 
ened  with  bounds  on  f'(x)  or  higher  derivatives. 

Lemma  U.2.6  If  F(x)  is  differentiable  s  times,  F(x)  is  symmetri¬ 
cal,  and  f^(x)  has  limit  zero  for  x->+«  and  for  x-^-», 
n*0,l,2, . . ,s,  then 


|B(x|x0)|- 


(U.10) 


n*0,l,2, . .  ,s 

f^(x)  is  the  n^  derivative  of  f(x) 

(See  Appendix  II  for  definition  of  P  (x)»-the  nth  Bernoulli 

n  function) . 

Proof :  Let  x^  be  the  Spearman  estimator  for  a  finite  number 
of  dose  levels,  x^*xQ+id,  i=0,±l,±2,  ,,.,±k.  From  (2.5): 


x,  *x  -d 
k  o 


k  l 
2(Pi-/2) 
-k  1  * 


E 


(\K)"Vd  f  4> 

-K 


Then  by  the  Euler-KacLaurin  formula  (see  Appendix  II): 
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From  (2,6)  and  (  .l)  it  is  seen  that  E(x|xq)«  lira  ECXjJXq)  , 

First  consider  the  limit  of  the  middle  two  terms  within  the 
brackets,  i.e., 

2  ^F(xQ+kd)-l+F(xo-kd)j  . 

Since  F  is  a  distribution  function,  the  limit  of  this  expression 
as  k  becomes  large  is  zero. 

The  first  integral  in  (It. 11)  can  be  rewritten: 


(Since  F(y)  is  assumed  symmetric,  F(y)  will  be  an  odd 
function  with  respect  to  y*p.)  Consider  the  limit  of  this 
integral  as  k  becomes  large. .  The  length  of  the  interval  of 
integration  remains  constant.  The  value  of  the  integrand 
approaches  ^  •  Hence, 

ip+kd+(xQ-p) 

[FCy^/ld^  »  ^2 

L  l!  d 

p+kd-(xQ-p) 
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Substituting  the  limiting  values  obtained  thus  far  in  ([*.11), 
the  expression  for  E(x|xo)  becomes: 

)»  lim  E(x,  [x  )*|i+  lim  d 
k— 4°°  0  k  — >  ® 

Thus  the  bias  can  be  bounded  as  follows: 


P1(x)f(xQ+xd)dx 


Before  taking  limits  the  integral  involving  and  f  could 
be  integrated  by  parts,  making  use  of  the  relationship: 


P^  (x)  -  (-1 )n“1  Pn-1(x; 

Repeated  integration  by  parts  would  lead  to  the  general 
expression  of  (U.10).  Q.E.D. 

U.2.7.  Expressions  for  the  Bound  Involving  Derivatives  of  the 
Tolerance  Distribution 

When  the  bound  (ii.10)  is  evaluated  for  n*0  and  n=l,  the 
results  are  respectively: 

|b(x)|  i  d/2 

\m)\  -  i- 

The  first  bound  is  identical  with  the  one  obtained  in 
paragraph  U.2,2.  without  the  assumption  of  symmetry.  The 
second  bound  is  of  the  same  order  in  d  as  that  obtained  in 
paragraph  U.2.U  but  the  constant  '/g  is  greater  than  /g. 
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Lemma  U.2.6  provides  a  sequence  of  bounds  on  the  bias 
in  increasing  powers  of  d.  For  example,  if  n  is  taken  to 
be  2,  the  expression  for  the  bound  becomes: 


lB<*M 


sup 

X 


®  sin  k  2TJ  x 
2 
1 


2^TT  3  k3 


d3  j  j  f"(x)| 


dx 


From  Appendix  II: 

*>  sin  k2fF  x 


sup 

X 


l 


•  .0080.  .  . 


A  simple  expression  for  the  integral  can  be  obtained  if  the 
density  f(x)  is  assumed  to  have  exactly  two  points  of  inflec¬ 
tion,  x*u±c,  with  the  absolute  value  of  the  derivative  of 

f(x)  at  these  two  points  being  f *  .  Then 

c 

p-c  i+C  ,<*> 

|f"(x)|dx  *\  f'*(x)dx  -  \  f"(x)dx  +  f"(x)dx 

p-c  p+c 

=  Uf' 
c 

Thus  the  bound  for  the  bias,  using  n*2,  becomes: 


Jb(x|xq)|S  *03 2d' 


(U.12) 


For  a  symmetrical  density  function  with  two 

points  of  inflection#  maximum  slope  f5  end  maximum 

c 

ordinate  f  : 

m 

f  -  In 

m  *  c 


Thus,  for  this  class  of  densities,  the  bounds  given  in  (luh) 
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and  (4.12)  can  be  combined  to  gL 
Theorem  4.2.7 

|B(i|x0)|<Kln.(.125d^  ,  .032d3fp  (4.12a) 

4.2.8  Computation  of  the  Bounds  on  the  Biaa  for  Some  Tolerance 
Distributions 

If  d  is  expressed  in  units  of  the  standard  deviation  or 
some  interpercentile  difference  of  the  tolerance  distribution, 
the  magnitudes  of  the  various  bounds  on  the  bias  obtained 
above  can  be  more  readily  compared.  Bound  (4.3)  becomes: 

|B(i|*0)|*  |  0  (p*f)  . 

Thus  to  assure  that  the  bias  is  less  than,  say,  10  percent 
of  the  standard  deviation  of  the  tolerance  distribution, 
d  should  be  less  than  20  percent  of  c.  ' 

If  the  bound  (4.4)  is  to  be  used,  the  modal  ordinate 
must  be  specified.  In  Table  4.1  bounds  on  the  bias,  computed  from 
(4.4),  are  given  both  in  terms  of  o  and  R.  R  is  the  distance 
from  the  20  to  the  60  percentile.  The  bounds  on  the  blew 
in  terms  of  a  and  R  for  the  various  distributions  are  quite 
similar.  For  the  four  distributions,  other  than  the  special 
one-particle  function,  a  choice  of  <K1.3o  will  assure  a  bias 
of  less  than  10  percent  of  o. 

The  bound  given  in  (4.12)  can  be  applied  to  the  first 
three  of  the  tolerance  distributions  listed  in  Table  4.1. 


Table  h.2  gives  the  bounds  for  the  three  tolerance  distri¬ 
butions  in  terms  of  both  a  and  R.  Note  from  Table  U.2  that 
if  d—  2,8.d  the  bias  of  the  Spearman  estimator  will  be  less 
than  10  percent  of  a  (using  i*.12). 

U*3  Effect  of  Dose  Mesh  Location  on  the  Variance  of  the  Spearman 
Estimator 

U.3.1.  Bound  on  jV(xjx0)-Vf 

Let  x  be  the  median  of  F(x)  (or  a  median  if  F(x)»i 
oO  “ 


does  not  uniquely  determine  x  .50)#  Then  F(tJ  (  l-F(t))  is 

non-decreasing  for  x-x  ^  and  is  non-increasing  for  x^x 

<  / 

Number  the  dose  intervals  so  that  xq*x  ^  x^.  Then 


L 

F(t)  |l-F(t)J  dt  -  dFt  ^1-Fj.J  i-0,-1,-2 


x.  , 
i-l 


i+1 

F(t)  { l-F(t)^  dt  -  dF,  f  1-F, 


i-l,  2,., 


F(t)  (l-F(t)j 


dt  «  d*.|*| 


Combining  these  inequalities: 


d2  F(xi)  [l-F(xi)l  +d/^  \  F(t)[l-F(t)Jdt 
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Table  li.l 


Bounds  for  the  Bias  of  the  Spearman  Estimator 
for  the  Infinite  Experiment  ^ 


y 


Tolerance 

Bound  for 

the  Bias 

Distribution 

as  a  Proportion 

of  a 

y 

as  a  Proportion  of  R 

.0567  ^ 

2 

u 

.0866^ 

>’2 

1 

3/ 

Logistic 

)  a 

R 

Normal 

.0U99  ^ 

>2c 

.0839 1 

d'2 

R 

Angular 

.0ii2?  ^ 

>2a 

.080I1  ^ 

d'2 

\ 

R 

Uniform 

.0361  ^ 

.0750^ 

f* 

R 

One -Particle 

.1250  r 

)  2c 

•1732 c 

y2 

R 

using  (U.U)  for  the  bound,  |  3(5c| xq )|  ^  f^ 


'  \ 

y 

2/ 

-  a  is  the  standard  deviation  of  the  tolerance  distribution 
and  d»  0  or 


T 


^  R  is  the  distance  between  the  20th  and  the  80th  percentile 


and  d-  p  ’  R. 


? 
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Table  U.2 

Bounds  for  the  Bias  of  the  Spearman  Estimator  for  the 
Infinite  Experiment 


Tolerance 

Bound  on  the  Bias 

y 

Distribution 

As  a  Proportion  of  As  a 

3/ 

Proportion  of  R  — ' 

Logistic 

.0101  p3cr 

.0568  j 

«  3 

R 

3/ 

Normal 

.0077  f3a 

•0529  | 

1 3 

P  R 

Angular 

.0075  P3u 

.0639  1 

?'3  » 

^ using  (U - 12 )  for  the  bound,  | B(x  |  x  ) j  *  .032d3f« 

o  tn 

is  the  standard  deviation  of  the  tolerance  distribution 
and  d*pa 

3  f  "bh  bht 

—'R  is  the  distance  between  the  2(j  and  80w  percentile  and 

R. 

Note:  Finney  (15,  16)  gives  the  actual  maximum  bias  for  the  normal 
and  logistic  distributions  for  various  values  of  d.Far-p  .of 
2-  B*.005a  for  the  normal  distribution  while  the  table 


gives  .062 a  as  the  bound;  at  p*3  Finney  has  B».107  compared 
with  a  bound  of  .208a  in  the  table. 
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Therefore  the  variance  (3.2)  satisfies  the  inequality: 

V(i|Xo)  -  Z  F(Xl)  [l-F(x.  )J  -  |  ^  F(t)[l-F(t)]dt- 

i  cm 

Similarly  an  upper  bound  for  the  variance  of  the  estimator 
can  be  obtained: 


V(x|xQ)  =  £  J  F(t)  jV-F(t)j  dt4^ 


Let  V  be  defined  as: 


V  -  |  J  F(t )  £  l-F(t) 


dt 


(U.13) 


Then 


li.3.2.  A  Tolerance  Distribution  for  Which  |v(x|x  )«*v|  Approximates 
the  Bound  Arbitrarily  Closely 

iu5) 

The  boundwDn  the  deviation  of  V(xjxo)  from  V  is  the 
Buprenum.  Consider  the  two  point  distribution  with  masses 
of  §  at  the  points  0  and  1.  Then  F(x)  has  the  form: 

F(x)  =0  x  <0 

F(x)  *  Vo  0  *  x<  1 


F(x) 


1  =  x 


Then 


V.|  \  F(t)  [l-F(t)| 


dt  » 


EH 


d  1 

EH  ’  3 


If  d*l-e,  1>  e>0,  then  it  is  possible  for  two  doses  to  lie 
between  0  and  1.  In  this  case. 
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v(Slxo)  ■  t~  (r?  *  *  §; 

•v(Slxo)-?l*  E,2_Erl 

v<Jlxo>-7  ■  e  (2t^) 

Therefore,  as  e  goes  to  zero  the  deviation  of  V(xjx  )  from 

p  I  o 

rf  a 

V  can  be  arbitrarily  close  to  +  ^  . 

Also,  using  the  same  example,  but  choosing  d  =  1+e, 


V(J|xQ)-?  -  -jf  •  ^ 


if  no  dose  levels  lie  between  0  and  1,  so  that  the  deviation 

d2 

can  be  arbitrarily  close  to  ~  • 

U.3-3.  Bounds  on  |v(:~:|x  )-V  |  in  Terms  of  Derivatives 


The  Euler-KacLaurin  f  uviulae  yield  better  bounds  on  the 
fluctuation  of  V(xjxQ)  due  to  the  placement  of  xc,  when  more 
stringent  conditions  are  imposed  on  F(x).  If  F(x)  has  a 
continuous  density  (and  the  first  moment  of  F  exists,  as  has 
already  been  assumed)  then  (see  Appendix  II): 

~  2  F(xQ+id)  [l-F(xo+id)j*  ~  j  F(xo+xd)  Jl-F^-ttd  )j  dx 


-  ~  \  P!(x)f(x0+xa)  Jl-2F(xo+*l)j 


V(x|xq)-?  =[^- j  P1(x)f(xo+xd)^l-2F(xQ+xd)Jdx[  (U.15) 


Integrating  by  parts  in  (U.15),  which  is  justified  if  F*  is 
continuous  and  if  the  resulting  integral  exists: 


V(xJxo)-V 


33  f(x0+xd)  (l-2F(x0+*d)) 


First  assume  that  F(l-F)  has  exactly  two  points  of  in¬ 
flection,  say  at  x=*c^  and  at  x*^.  Then  the  integrand  is 
positive  for  x  <"c^  and  for  x  )  and  negative  for  x  between 
and  c^.  Then 

If  F  is  symmetrical, 

fv(x|xo)-7|  4  ^  f(Cl)  [l-2F(cx)] 

Also 

lv(£lxo)-r|  ’  35  (‘‘.W) 

li . 3 - U  Computation  of  Bounds  on  jv(x|xo)-v|  for  Some  Tolerance 
Distributions 

To  illustrate  the  magnitude  of  the  bounds  on  |v(xjxo)-7( 
express  d  in  units  of  the  standard  deviation  or  an  inter¬ 
percentile  deviation  of  F(x)  and  express  j V(x|xo)-tf|as  a 
proportion  of  V.  Thus,  using  (Lull;)  and  (U .16 )  ; 


lvg|*0)-d  \ 

V  j  F(x)[l  -F(x)j  dx 

P<x|xo)-7|  §2fm 


(Lul6a) 
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See  Table  U.3  for  numerical  examples, 
if  d  is  less  than  .Ua  then 

V 


From  the  second  column. 


<  nr.  ,  for  each 

*  20  percent 


of  the  five  tolerance  distributions.  The  fourth  column  indi¬ 
cates  that  for  the  logistic,  normal  and  angular  distributions, 

a  d  less  than  .8 cr  assures  that  |^(XIXQ)  ^  <  _ 

-  »  20  percent. 

V 

5.  .RANDOM  LOCATION  OF  THE  DOSE  MESH  IN  THE  INFINITE  EXPERIMENT 
5.1  Introduction 

Random  location  of  the  dose  mesh  is  accomplished  for  the 
infinite  experiment  by  fixing  the  dose  interval  d  and  randomly 
choosing  the  dose  level,  x  ,  from  the  uniform  distribution 
over  the  interval  (0,d).  Even  though  no  effort  is  made  to 
randomly  locate  the  dose  mesh,  in  certain  routine  screening 
procedures,  at  least,  the  tolerance  distributions  are  essen¬ 
tially  randomly  located  with  respect  to  the  fixed  dose  mesh. 


5.2  Unbiasedness  of  the  Spearman  Estimator 

The  expected  value  and  variance  of  the  Spearman  estimator 
for  random  choice  of  xq  will  be  denoted  by  E(x)  and  V(x). 

Irwin  (22)  and  Finney  (15,  16)  pointed  out  t  hat  when  the 
location  of  the  dose  mesh  is  selected  at  random  the  Spearman 
estimator  is  unbiased.  This  is  shown  as  follows: 
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Table  U.3 


Bounds  for  the  Relative  Deviation  of  V(x|xq)  from  7 

|V(5K)  -  7| 


Tolerance 

Distribution 


Bound  on 


First  Bound  — ' 
“2/ 


Second  Bound 


TT 


IT 


- — ^ - Tr 

.27U1  D2  or  .6Uoh 


Logistic 

.U535p" 

» 

or 

.6932 

Normal 

.UU31-  ^ 

or 

.7U58p 

Angular 

-U3U7  P 
\ 

or 

\ 

.81830 

Uniform 

•  U330  p 

t 

.5000  ^ 

or 

\ 

.9000 p 

One-Particle 

or 

\ 

.6930^ 

V 

.2357^)' 

.1981 


or  .6677 


or  .7020 
Not  applicable 
Not  applicable 


y 


y 


3/ 


The  first  bound  is  computed  from  (U.lUa)  and  the  second 
bound  from  (U.l6a).  (See  Table  5*1  for  the  values  of  7 
for  the  several  distributions.) 

is  the  ratic  of  the  <~ve  interval  to  the  standard  devia¬ 
tion,  d=^)  a. 

p '  is  the  ratio  of  the  dose  interval  to  distance  (R) 
from  the  20^  to  'the  80tn  percentile,  d-D'I  . 


‘f 


Q —  cu~ 
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rd 

E(x)  =  i  E(x|xq)  dxc 


(5.1) 


l 

d 


d 

00 

■  X  +id+d 

2  (x  +d/  +id) 

0  c 

— <x> 

O 

|  dF(x) 

0 

J 

'  x  +id 

dx 


kz 

— oo 


co 

J* 

—  CO 


(x  +id+d/p) 

0  c- 


x  +(i+l)d 
o 

dP(x) 


0  l 

J  . . 
x  +id 

0 

f 'i*l ;d+d/. 

• 

u+<V 2 

,  u 

dF(x) 

- 

d  / 

u-  A 

dx 


du 


1 

d 


1 

d 


^+d/2 

u  |  dF(x)  du 


[u"  '2 
t  dF(j 

J  d/ 

J  u-  / 


x+d/2 

udu  dF(x) 
d  / 

x-  4 


xdF(x)  =  p 


5 . 3  Ihe  Mean  Square  Error  of  the  Spearman  Estimator 

Let  E  denote  the  expectation  with  respect  to  x  over 
o 

interval  (0,d).  Let  B  (jcjx  )  denote  the  bias  of  x  given  xq. 
The  mean  square  error  of  the  estimat  or  is  the  variance  and 
can  be  written: 


V(x) 


(x-p)2[xo 


V(xj 


Xo> 


B1'  (x 
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E 


|v(*|xo)]  *EX  [b2(xxo)|  15.2) 


The  first  component  of  V(x)  in  (5.2)  can  be  evaluated: 


Ex  V(xjxo) 

n  - 


2  » 

i-  s  F(xo+id)  [l-F(xo+id)] 


d 

F(xQ+id)  j^l-F(xo+id)  j  dxQ 
o 

xQ+(i+l)d 

F(x)  |  1-F(x)l  dx 

x  +id 
o 


dx 

o 


*= 


d 

n 


F(x)  [l-F(x)j  dx 


(5.3) 


Note  that  (5.3)  is  the  same  as  (U.13),  denoted  by  7.  Thus 
V  is  the  average  of  the  conditional  variance  of  x,  taken  over 
the  location  of  the  dose  mesh, 

A  simple  expression  for  the  second  component  of  V(x)  has 
not  been  obtained.  The  variance,  V(x),  is  written: 


V(X)  =  a 


(x)  |  1-F(x)j  dx+E^ 


B2(i 


(S.U) 


Using  the  bound  for  the  bias  over  values  of  xq  in  ( !*•  3 )  “ 

V(x)  =  V  +  0(d2) 

For  tolerance  distributions  satisfying  the  conditions  given  in 
section  U.2.7: 


V(x)  =  V'  +  0(d6) 


The  second  component  of  V(x)  (5.U)  is  independent  of  n  and 
of  smaller  order  in  d  than  the  first  component,  ?  .  V  con¬ 
tains  d  only  in  the  form  of  theihctor  d/n>  i.e.  the  inverse 
of  the  number  of  subjects  tested  per  unit  interval  on  the 
dose  scale, 

5. U  Values  of  V  for  Severs!  Tolerance  Distributions 

An  approximation  to  the  variance  of  the  Spearman  esti¬ 
mator  for  the  case  of  a  normal  tolerance  distribution,  given 
by  Gaddum  (19),  is  equivalent  to  V  ,  Finney  (15,  16)  com¬ 
puted  V  for  the  normal  and  logistic  distributions.  Table  5.1 
' gives  the  values  of  7,  as  proportions  of  a  and  as  proportions 
of  the  distance  from  the  20^  to  80^  percentile  (R),  for 
several  tolerance  distributions  (see  Table  2.2  for  defini¬ 
tions  of  these  tolerance  distributions.) 

In  section  U.3  it  was  seen  that  V  can  be  regarded  as  a 
good  approximation  to  V(x|:„o>,  and  that  V(x|xq)  can  deviate 
from  V  by  at  most  . 

In  the  present  part  V  was  seen  to  be  less  than  the  uncondi¬ 
tional  variance  of  x  ,  where  the  error  is  slight  if  d  is 
small.  In  the  next  section  V  will  be  established  as  the 
asymptotic  variance  of  x,  as  defined  in  the  same  section. 

6.  LARGE  SAMPLE  PROPERTIES  OF  THE  SFEAHM4N ESTIMATOR  FOR  THE 
INFINITE  EXPERT!  ^NT 

6.1.  Large  Sample  Experiments 

The  experimental  design  for  the  infinite  experiment 
consists  of  fixim.:  two  numbers:  the  number  (n)  of  subjects 
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Table  5.1 


Values  of  V 

for  Several 

Tolerance 

Distributions 

Tolerance 

,y 

Distribution 

Logistic 

.5513^ 

2  y 

a 

n 

1.2831  £ 

2 

n 

Normal 

2 

>£- 

n 

1.5983  j3 

R2 

»  '  £- 

n 

Angular 

.575od 

2 

c 

rT 

2.0376  j! 

>.  5! 

n 

Uniform 

.^7?Up 

2 

a 

n 

2.h9h2 

si 

n 

One -Particle 

.5000^ 

£_ 

n 

.9601^ 

R2 

>«  — 

n 

y  d  i 

V  -  - 

n  ' 

F(l-F )dx  ; 

these  values  were  used  in  Table  ii.3 

■ 

mCO 

2 / 

P  is  the  ratio  of  tne  dose  interval  to  the 

standard 

deviation,  d*flo 


3/  th 

p  '  is  the  ratio  of  the  dis+ance  (H)  from  the  20  to  the 

80^  percentile ,  d*  ^  '  R 
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to  be  tested  at  each  of  the  levels,  and  the  distance  (d) 
between  dose  levels.  Denote  such  a  desfn  by  D(n,  d). 

The  large  sample  experiment  Is  usually  described  in  terms  of 
a  fixed  d  and  increasing  n.  For  the  Spearman  estimator  this 
method  of  increasing  the  else  of  the  experiment  will  not 
yield  a conoi stent  estimator.  The  second  component  of 
v(x)  in  (5.4)  involves  the  conditional  bias  independent  of 
n,  so  that  V(x)  does  not  go  to  zero  as  n  goes  to  infinity. 
This  points  up  the  need  for  a  more  general  concept  of  large 
■ample  experiments. 

Let  n'  denote  the  average  number  of  subjects  tested  per 
unit  on  the  dose  scale,  (n?»n/d) .  Then  the  large  sample 
experiment  is  obtained  by  letting  n*  go  to  infinity.  The 
choice  of  the  corresponding  voJues  for  n  and  4  will  be  made  to 
minimize  the  mean  square  error  for  fixed  n* , 

6.2  Optimum  Choice  of  n  and  d  for  the  Spearman  Estimator 

For  fixed  nf  and  random  choice  of  xQ  the  variance 
of  the  Spearman  estimator  is  minimized  by  choosing  n  and  d  as 
small  as  possible,  i.e,,  n»l  and  d*Vni»  This  follows 
from  the  following  theorem: 

Theorem  6.2.1:  The  mean  square  error  of  the  Spearman  estlmater 
based  on  groups  of  n  subjects  tected  at  dose  levels  d 
units  apart  is  greater  than  the  mean  square  error  for 
single  subjects  tested  at  dose  levels  */  units  apart. 

Proof:  Denote  the  two  mean  square  errors  by  MBEQ  and 

respectively.  Denote  the  corresponding  biases  conditional 
on  xQ  by  \(Xq)  end  b^x^)  respectively. 
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*®n  -  5  /'(1-F)<U  +  E  t^(x0)) 

0 

(6.1) 

**1  ■  ^rfr(i-T)ix  +  E  tt>|(x0)  1 

(6.2) 

The  first  terms  on  the  right  hand  sides  of  (6.1)  and 
(6.2)  are  Identical.  Therefore  it  must  he  shown  that 


yw>  >  ywi 


(6.3) 


The  conditional  biases  can  be  written: 


00  rU 

t>_(xj  *•  £  /  (x0+id+d/2-x)dF 

i—co  J 


xQ+id+d 


Vid 


.  x  +id+^i^: 
od  n-1  pO  n 

V‘“f 

The  left  hand  side  of  (6.3)  is:; 


(6.4) 


(6.5) 


Vbn<*0» 


V4 


l/  ^*0 


+  (e+l)d 

3  .So  J  Vxo)dxc 


n-1  rO 


n-1 


*0*  £ 

vf 


3  /  bn(*0-£  >*b 


(6.6) 
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The  expression  b  (Sq+  |j=)  appearing  in  the  Integrand  of 
of  the  right  hand  side  of  (6.6)  can  he  rewritten  in 
terms  of  b^(x^) .  From  (6.U): 


x.  +  +  id  +  d 

u  f  n 


*>n(x0  + 


oo 

E  /(*0  +  T  +  id  +  d/2  -  x)dF 
i«  -c>  * 


(6.7) 


x+  ~  +  id 
o  n 


S  «  0,1, 2, « • • (n-1) « 


oo  n-1  r  . 

E  E  /  (jl  +  |p  +  id  ♦  d/2  -  x)dF  (6.8) 
i-  -Oj  J»0j 


* 0  +  a  + « +  oate 


*0  ♦  T + 14  +  ir 


Then,  adding  and  subtracting  in  the  integrand 

of  (6.8)  : 


*o  ♦  id  +  Usss^l 


V*0 +  r>  -  ?  “£/[*< 0 + 11  ♦  5-  - xldr 

i*  -co  J«0  w 


xn+id+ 
u  n 


d  oo  n-1 

+  I  iS  -oo  j£o 


Xo  + 14  +  Lxzrn 

J  (i  - 

,  +  id  +  Llaah 


(6.9) 


From  (6.5),  expression  (6.9)  can  he  written: 
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n-1 


VV  ¥  >  *  VV  T>  +l  A'1'  ^Aj„  <6-10> 


J*o 

xo  +  id  +  0±S±lU 

00  r  n 

where  A.  _  ■  E  /  dF 

J+<  -<0  ♦  14  ♦ 

u  n 


n-1 

Note  that  E  A 

*■0 


1  for  any  J 


(6.11) 


Since  b^  is  periodic  with 

(6.10)  that 

V  V  T  ) 

*  W  ♦ 1 

n-1 

E 

J=0 


B+J 


(6.12) 


Substituting  (6.12)  in  the  right  hand  side  of  (6.6): 


t*n(V> 


d 

2 


n-1 

E 

>0 


(6.13) 


On  squaring  the  expression  in  the  integrand  in  (6.13)  the 
middle  tern  will  be  : 


(6.14) 


n-l 

£ 

S-0 


J%  (1-  ^ 


Summing  fir  at  on  e  the  sum  of  the  Is  one  and  the 
sum  with  respect  to  J  will  then  he  zero.  Therefore 
(6.13)  becomes: 


Vbu( V  1  ■  B*  •bI(xo)1  + 

0  \) 


*  I 


Pi- 


PH-1 

/  E 

J  mmO 


n-l 

(  t  (1 

a*0  j«0 


(6.15) 


The  second  tern  on  the  right  hud  side  of  (6.1?)  csnnot 
he  zero  If  F  is  a  distribution  function.  Therefore, 


hi 

6.3  Large  Sample  Properties  of  the  Spearman  Estimator 

If  n'  ia  increased,  with  n»l  and  d*Vni>  «• 
required,  it  follove  from  (5.^)  and  the  fact  that  the 

p 

bias  squared  has  a  bound  of  order  d  ,  that  the  estimator 
is  consistent  and  the  variance  is 

V(i)  -  |,y'F(l-P)ax-tO(1/n,2)  (6.9) 

-00 

It  is  convenient,  in  the  case  of  a  sequence  of 
random  variables,  to  approximate  the  variances  by 
simpler  terms  correct  to  order  n"  .  (When  the  sequence 
does  not  have  variances,  the  variances  of  a  sequence  of 
limiting  distributions  may  be  used.)  Such  a  sequence 
of  approximations  will  be  called  the  asymptotic  variances. 
In  this  sense  the  first  term  on  the  right  of  (6.9)  will 
be  called  the  asymptotic  variance  of  x  and  will  be  denoted 
*y  vA(x). 

vA(i)  - 

-co 


(6.10) 
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7.  URGE  SAMPLE  EFFICIENCY  OF  THE  SPEARMAN  ESTIMATOR  FOR  THE 
INFINITE  EXPERIMENT 

7.1  Previous  Comparisons  of  the  Spearman  Estimator  with  the 
Maximum  Likelihood  Estimator 

Finney  (l£,l6)  computed  the  asymptotic  variances  of  the 
maximum  likelihood  estimator,  averaged  over  choices  of  xq, 
for  the  finite  experiment,  for  the  normal  and  logistic  tol¬ 
erance  distributions,  and  then  took  the  limit  as  the  number 
of  levels  went  to  infinity.  He  compared  these  values  with  the 
oean  square  error  of  the  Spearman  estimator  over  choices 

of  x  for  the  same  two  distributions.  The  ratios  were  .9811; 
o 

and  1.0000  respectively. 

Cornfield  and  Mantel  (10)  showed  that  for  the  logistic 
tolerance  distribution,  the  maximum  likelihood  estimator 
and  the  Spearman  estimator  were  approximately  equal  and  this 
algebraic  approximation  improved  as  d — >0.  Bross  (9)  evalu¬ 
ated  some  sampling  distributions  through  enumeration  for  the 
maximum  likelihood  estimator  ana  the  Spearman  estimator.  He 
used  the  logistic  tolerance  distribution,  four  dose  levels, 
with  n*2  and  also  n*{>.  In  «.l  1  cases  examined, 
the  Spearman  estimator  was  concentrated  mors  closely  about 
the  true  mean  tolerance  than  was  the  maximum  likelihood 
estimator.  These  computational  results  were  reproduced  by 
Haley  (21)  for  the  normal  tolerance  distribution. 


These  results  indicate  that  the  Spearman  estimator 
compares  favorably  in  precision  with  the  maximum  likelihood 
estimator,  at  least  for  the  normal  and  logistic  distributions 


In  this  section  the  asymptotic  efficiency  of  the  Spearman 
estimator  is  defined  and  various  results  are  reported  con¬ 
cerning  tolerance  distributions  that  minimize  or  maximize 
this  efficiency.  Efficiencies  for  the  common  tolerance  dis¬ 
tributions  are  given,  the  v  alues  for  the  logistic  and  normal 
being  the  same  as  the  ratios  given  by  Finney. 

7.2  Definition  of  Asymptotic  Efficiency 

The  efficiency  of  an  estimator  con  be  defined  in 
terms  of  the  quantity,  I,  called  the  information: 

i  -  e(-^~-£)2 

o  *-*• 


where  f  is  the  frequency  f ’unction  for  the  random  variables 
on  which  the  estimator  depends  (11). 

For  the  infinite  experiment  with  random  choice  of  xq, 
the  information  is  (see  Appendix  I): 
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n 

d 


— oo 


FU7[  X-FlxJJ^ 


The  asymptotic  efficiency,  E,  of  an  estimator  for  the 
infinite  experiment  will  be  defined  as  the  r  atio  of  to 
the  asymptotic  variance  of  the  estimator.  For  the  Spearman 
estimator. 


hh 


V, 


VA(x) 


F(x)[l-F(x)]  dx 


F  2(x) 

-J i_ 


J  F(x)[l-F(x)] 


dx 


Vj  is  the  asymptotic  variance  (n— >■»)  of  the  maximum 
likelihood  estimator,  so  that  E  measures  the  efficiency  of 
the  Spearman  estimator  relative  to  the  maximum  likelihood 
estimator. 

7.3  The  Spearman  Efficiency  for  Several  Tolerance  Distributions 
The  f  ollowing  sections  present  computational  results 
for  specific  tolerance  distributions.  The  results  are 
summarized  in  Table  7.1. 

7.3.1.  Logistic 

The  efficiency  for  the  logistic  is  1.0  since  the  logis¬ 
tic  distribution  satisfies 


P  F(l-F)  -  Fx 


7.3.2.  Normal 

Finney  evaluated  the  efficiency  for  the  normal  tolerance 

F  2 

distribution.  Th.  integral  of  ^  t0  ^  obtained 

through  numerical  methods.  The  result  is  E*.98lli 
7-3.3.  Angular 

For  the  angiilar  distribution,  both  integrals  are  easily- 


evaluated  and  £*.8106 


7.3. U.  Uniform 


U5 


The  definition  of  E  is  not  applicable  to  the  uniform 
distribution  since  the  regularity  conditions  used  in  obtain* 
ing  I  are  not  fulfilled.  The  particular  difficulty  is  that 
the  distribution  function  is  not  differentiable  for  all  values 
of  p  (see  Appendix  I). 

7.3*5.  One  Particle 

The  integrals  for  the  one-particle  distribution  are; 


d 
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F(l-F)dx  -  ^  |  a 


n 
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F(l-F) 


n 

dx-  -  a 
d 


-y 

y^e  dy  ^  n>a“  1_  _  2,l40U3na 
^‘y  d  lk3  d 


2 

r  , 

l-e 


o''  o 

and  the  efficiency  is  .83191... 

7.3.6.  Algebraic 

The  above  examples  all  involve  distributions  with  first 
moments  for  all  values  of  the  parameter.  The  algebraic 
distribution  (Table  2.2)  does  not  have  a  first  moment  for  all 
values  of  the  one  parameter: 


F(x;  s)  *  l-x’S  s>l,  x  =  l 

s 

This  distribution  has  a  first  moment  with  value  u  *  — r 

^  s-1 

if  s>  1.  For  s  =1  the  first  moment  does  not  exist. 


d 

n 


F (l-F)dx  -  - - - - 

n  (s-l)(2s-l) 


s>  1 


1 
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n 

d 


f 


—40 


F 

— ti -  dx  * 

F(l-F) 


J  2(5-1)^  2  - 

k=l  (ks-1)3 


s>  1 


and 


E  =• 


2-1/, 


~T“ 


k=li  k-  /s 


Thus,  in  contrast  with  the  preceding  examples.,  E  depends  on 


s. 


lim  E  =  i  lim  (2-Vs) 
s — s — *1 


1 


lim  2 
s — *L  k=l 


(l-X/s)3 

(K-Vs)3 


The  lirr.  -  in  +he  denominator  is  one  since: 


_  00  (s-l)3  *  (s-1)3 

lim  2  -i -SZLL »  lim  (1  +  2  -  )  =*  1 

s— >1  k=l  (ks-1)3  s— *1  k=2  (ks-1)3 

Also 

lim  E  -  *=  .8319... 

s  — *°°  1 


Note  that  the  limit  of  the  efficiency  as  s  becomes 
large  is  identical  with  the  efficiency  for  the  one-particle 
tolerance  distribution.  This  might  be  anticipated  by  re¬ 
writing  the  two  distributions  from  Table  2.2. 

Algebraic  Distribution:  l-(x+l)“s  x  =  0 

One-Particle  Distribution:  l-e'"^s'’3^x  x  «  0 


1 

^ — J 


i 

“’in 


1*7 

The  ratio  of  the  i^*1  moments  of  the  two  distributions  is 
(s-l)(s-2).. .(s-i) 

-TT-  -  '  ■  r  ■  -  i  "  , 

( s~l )i 

This  ratio  goes  to  one  as  s  becomes  large. 
l.h  The  Logistic  Tolerance  Distribution  and  the  Spearman  Estimator 
Theorem  7.U.1: 

The  logistic  tolerance  distribution  is  the  only  symmetri¬ 
cal  tolerance  distribution,  with  a  translation  parameter  as 
the  single  'unknown  parameter,  for  which  the  Spearman  estima¬ 
tor  has  full  efficiency. 

Proof;  Let  \x  be  the  translation  parameter  and  let  the  toler¬ 
ance  distributions  be  written  F(y-p).  The  efficiency  for  a 
given  distribution,  F,  is: 


Let  G  be  an  extremal  function  (symmetrical,  differentiable) 
of  the  functional  Y(F)*  j^E(F)j  \  Let  V(x)  be  any  function 
satisfying  the  conditions: 

V(x)*V(-x)  (7.1) 

V(x)  is  differentiable  for  all  x  (7.2) 

for  all  t  in  a  neighborhood  of  t*0,  G(x)+tV(x)  is  a 

distribution  function  with  first  moment.  (7.3) 
Then  y(t)=Y  j^G(x)+tV(x) 
at  t*0,  and  y’(0)»0. 


is  a  function  of  t  differentiable 
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Table  7.1 


Large  Seaple  Efficiency  (E)  of 

the  Spearman  Estimator  for  the 

Infinite  Experiment  for  Several 

Tolerance  Distributions  ^ 

TOLERANCE  distribution 

EFFICIENCY  -/ 

1.  Logistic 

1.0000 

2.  Normal 

.9814 

3*  Angular 

.8106 

4.  One-Particle 

.8319 

5*  Algebraic 

.500  *  E  ^  .8319 

2/  See  sections  7.3*1/  7*3*2,  7*3*3,  7*3*5  and  7*3*6  for 
computations. 


-  "ti—)  -  -.)>■  1  -ag 


U9 


Let  G’(x)  be  denoted  by  g(x)  and  V’(x)  be  denoted  by 
v(x).  Then 


y'  (o)* 


2 

__s _ 

G(l-G) 


dx 


V(1-2G)  dx 


g2(i-g)2 


G(l-G)  dx 


J 

Since  G  has  a  symmetrical  density  and  V  is  symmetrical,  the 
integrand  of  the  following  integral  is  an  odd  function  and 


2G(1-G)gy 

G2(l-G)2 


dx 


0 


Then 


y»(o)«  1  V(1-2G) 


c,« 


cri? 


G  (1-G) 


dx 


where 


G 


1" 


j  G(l-G) 


—  CD 


G(l-G)dx 


V(x)  can  be  any  function  satisfying  conditions  (7.1),  (7.2) 
and  (7.3)  and  l-2G(x)  cannot  be  identically  zero  on  the 
infinite  interval.  Therefore  the  necessary  condition  for 
G(x)  to  be  an  extremal  function  of  Y(F)  is  that 


5o 


'T 


02g2(x) 


=  0 


G(m)  |l-G(x)]|2 
This  implies  that  G(x)  is  of  the  logistic  form: 


G(x)  - 

Q . E.D • 


7.5  Distributions  vith  Efficiency  of  the  Spearman  Estimator  Close 
to  Zero 

If  F(x)  has  a  first  moment  the  variance  of  the  Spearman 
estimator  exists,  i.e.  the  integral  with  respect  to  F(l-F) 
is  finite.  Then,  if  the  information  is  finite,  the  efficiency 
of  the  Spearman  estimator  is  greater  than  zero. 

Distributions  with  a  translation  parameter  as  the  single 
unknown  parameter  can  be  specified  for  which  the  Spearman 
efficiency  is  arbitrarily  close  to  zero.  Consider, 

f  x 


Fe(xi  n)  *  K(e) 


dt 


i  [- 


lipr 

l+2e 


1+e 


-<  x<  •> 
-><  - 
0<  e 


The  efficiency  can  be  made  arbitrarily  small  by  choosing  e 
close  to  zero;  V(e)  is  unbounded  as  e  goes  to  zero,  while 
1(e)  is  bounded  away  from  zero. 

Consider  the  f  oil-owing  bound  for  V(e): 


V(e) 


dx 
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-  |  2j  Fe(x;o)  j\-Fe{x;o)  j  dx 

o 

■  n  |  [l-Fe(x,-o)]dx 
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>  d 

'•  /* 

|  K(e)dt 

n  ! 

7  ?! 

l+e 

j  j  i* 

o  x  |  l+2e  J 

> 


dx 


> 


dX(e) 


1  x 


dt 


dx 


>  dK(e) 


1  1_ 
l+2e  2e 


(7.U) 


The  constant  K(e)  necessary  to  make  F  a  distribution  is 
greater  than  ±  for  e < 1;  therefore  V(e)  goes  to  infinity 
as  e  goes  to  zero. 

1(e)  can  be  written: 


1(e) 


* 

c 


r  I  S-Lle 

0  JJ. _ 


j  Fe(xilx)  |l-Fe(xj^)J 


dx 


To  show  that  1(e)  is  bounded  away  from  zero  as  e  goes  to 
zero,  first  note  that 


Using  the  symmetry  of  the  integrand  and  F(l-F) 

r  •% 

I  1^1  2 

1(e)  -  2 


<  3. 

al 


-LKJdi —  dx 
F(l-F) 


K2(e)dx 


1  [*  3b] 


5  1  5*5e 


Then  the  following  inequalities  are  obtained: 


1(e)  «  8K2(e) 


[  h ) 


dx 

2  1  2+& 


-  8K2(e)  -  1 


,2+2e 
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-  8K2(e)  — 

2^+^e  3+Le 
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(7.5) 


3  2  3+be 

It  follows  from  (7.1:)  and  (7.5)  that  E(e)  goes  to  zero  as 
e  goes  to  zero. 

7.7  Two  Parameter  Families  of  Tolerance  Distributions 


The  r  esults  of  the  previous  paragraphs  are  applicable 
without  modification  to  the  case  of  scale  parameter  unknown 
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when  the  tolerance  distribution  is  symmetrical.  (Estimation 
of  the  scale  parameter  itself  isduscussed  in  Appendix  III.) 

Let  the  tolerance  distribution  be  of  the  form  F(y) 
where  y~^(x-p).  Let  both^  and  p  be  unknown.  The  infinite 
experiment  information  matrix  is  given  in  Appendix  I. 

Letting 

l  ' 

tWdt 


*1 


Wdt 


(t-ATit 


Wdt 


where  W(t) 


[f*  (t)] 


F(t)  [l-F(t  j] 

,11 


The  inverse  element,  I  ,  corresponding  to  p  is 


d  a/*  A2 


r  cc 

„r 

j  W(t)dt 

Note  that  if  A^*0,  then  I^^/j  when  I  is  the  information 

for  scale  parameter  known.  If  F  is  a  symmetrical  distribution 

then  W(t)  is  $rmr.etrical  and  A^*0. 

Note  also  that  if  A-^/O,  the  effect  is  to  increase  the 
11 

value  of  I  above  the  corresponding  value  for  s  cale  parameter 
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knovn.  The  variance  of  the  Spearman  estimation  ia  unchanged. 
Hence,  the  efficiency  for  the  Spearman  estimator  would  be 
greater  in  such  cases  for  scale  parameter  unknown  than  for 
scale  parameter  known. 

8.  THE  SPEAHMAil  ESTIMTOR  FCR  THE  FINITE  E]<PERIMENT 


8-.1  Finite  Experiments 

In  previous  sections  (3>  4,  6,  7)  the  range  of  experi¬ 

mentation  was  infinite.  The  results  obtained  are  useful  in 
designing  and  interpreting  experiments  in  which  the  dose 
levels  cover  the  greater  portion  of  the  range  of  F(x),  say 
from  .01  to  .99.  As  a  supplement  to  these  results  it  is  of 
interest  to  investigate  the  effect  of  using  a  finite  set  of 
dose  levels . 

Lot  xq  be  an  a  priori  estimate  of  [i  and  let  the  experi¬ 
ment  involve  2k+l  dose  levels  regularly  spaced  over  the  in¬ 
terval  (x  -a,  x  +a).  The  dose  levels  arc  x.*x  +id,  i-0, 

±1,  ±2,...±k,  with  kd*a.  Let  N  be  the  total  number  of  sub¬ 
jects  used  in  the  assay,  N=(2k+l)n. 

8.2  The  Spearman  Estimator 

The  variance  and  bias  of  the  Spearman  estimator  are 


k 

V_(*|  xj.2-  Z  F  (1-F. )  (8.1) 


a  l  o  n  i'  x' 


a  k_1  rt 


Ba(iK>  ■  (va_  4)F(va*  4)+  l  cv  4Hri+i_Fi) 
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*(x0*a*  /2) 


1-F(x  +a-£) 
o  c 


-M- 


(8.2) 
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Theorem  6.2.1  Let  the  range  of  dose  levels  be  (xo-a,xo«*),  tho 
total  sample  size  be  N,  and  the  numbers  of  subjects  at  each 
dose  level  be  equal  (n).  Then  the  maximum  variance  of  x  over 
all  possible  F  is  minimized  by  minimizing  d  (i.e.  by  maximiz¬ 
ing  the  number  of  dose  levels). 

Proof:  From  (8.1)  the  variance  of  the  astimator  is 


V5 


,  ,a.2  (2k*l) 

V  V 
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Un 


The  bound  is  attained  for  F(x)  a  two  point  distribution 
defined  b> : 

P(x*x  -a-e)  *  Vo  e>  o 

o  '  2 

P(x=XQ+a+e)  *  Vg 

The  bound  is  minimized  by  choosing  k  as  large  as  possible, 
i.e.  by  choosing  d  as  small  as  possible.  y.E.D. 

There  are  distributions  for  which  an  increase  in  k  results 
in  an  increase  in  the  variance  of  x.  Consider  the  distribu¬ 
tion  given  by: 

p(x'v34a)'1/2 

p'-v34a)=1/2 


K 
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The  variances  for  three  and  five  levels  are: 


k*l  :  Va(x  j  xq) 


UN 


k-2 


V  (x  x  ) 
a '  o' 


15a: 

16N 


The  variance  for  the  normal  distribution  for  several 
values  of  k  is  shown  in  Table  8.1. a.  The  results  are  for 
xo*n.  For  the  dose  ranges  used,  the  variance  decreases  as  k 
increases . 

For  asymmetrical  location  (xq/pl)  the  Spearman  estimator 
will  be  biased.  Table  8.1.b  presents  the  mean  square  error 
for  the  normal  distribution  for  several  values  of  a,  k,  and 
N,  for  several  values  of  x  . 

8.3  Information  for  the  Finite  Experiment,  Scale  Parameter  Known 

Denote  the  information  for  the  finite  experiment  described 

in  section  8.1  by  I  (x  ).  Then 

a  o 


I  (x  )■ 
a  o 


2k +1  -k 


F 

p,  k _ 

f(t} 


Table  8. 2. a  presents  the  values  of  2T  /  \  for  the 

a  Ia(xo) 

normal  distribution,  for  x  *ji,  for  several  values  of  a  and 
k.  The  results  show  that  there  are  extreme  situations 
(three  levels  placed  at  q.-10o,  (i,  and  (i+10a)  for  which  an 
increase  in  k  results  in  a  decrease  in  the  information. 
However,  when  a  is  2a  or  less,  the  denser  the  dose  levels, 
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Table  8.1. a 

Variance  (V  )  of  the  Spearman  Estimator  for  the  Mean  of  a  Normal 

3 

Tolerance  Distribution,  for  the  Finite  Experiment 


Middle 

Dose 

location 

Dose 

Range 

NVft(x  |  x 

- 1 - 

o 

o> 

X  -LI 

(  °  ) 
v  a  ' 

k-1^/ 
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k»3 

k-li 
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k«20 

k-UO 

0 

1/2 

.51 

.36 

•  32 

•  30 

.27 

.25 

,2ii 

0 

1 

1.55 

1.18 

1.06 

1.00 

.91 

.87 

.83 

0 

2 

3-5U 

2.81 

2. $8 

2.51 

2.36 

2.25 

2.19 

0 

h 

12.00 

5.89 

5.08 

5.02 

U- 75 

U-  58 

a. 50 

0 

10 

75.00 

31.25 

19.  hh 

15.18 

12.37 

11.  U6 

11.23 

The  number  of  dose  levels  is  2k +1. 
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Table  8.1.b 


Mean  Square  Error  (MSE^)  of  the  Spearman  Estimator  for  the  Mean 
of  the  Normal  Tolerance  Distribution,  for  the  Finite  Experiment 


Middle 

Dose 

Location 

x  -p 

<~3-> 

Dose 

Range 
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M  [  KSEa(i  |xo) 
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- - 

k«l— ^  k*2  k*b 

0 

±1 

1.55  1.18 

1.55  1.18 

.5 

1.51  1.22 

1.92  2.36 

1.0 

1.55  1.63 

b.5l  8.60 

2.0 

5.13  7.99 

b7.09  77.35 

3.0 

23.30  31.20 

232.36  3H.65 

b.o 

1 

62.5?  75.67 

625.65  756.66 

0 

±2 

3.5b  2.81 

3.5b  2.81 

.5 

j 

3.38  2.79 

3-38  2.79 

1.0 

j 

3.22  2.71 

3.22  2.76 

2.0 

3.29  2.37 

3.b8  5.38 

3-0 

2.6b  5.b5 

11.86  b7.bl 

U.o 

| 

11.20  23.35 

109.60  232. bl 

5.o 

t 

I 

liO.  12  62.57' 

b01.06  625.66 

0 

j 

±b 

5.89  5.08 

5.89  5.08 

.5 

5.6b  5.08 

5.6b  5.08 

1.0 

5.39  5.08 

5.39  5.08 

2.0 

5.89  5.0J 

5.89  5.07 

3.0 

\ 

5.37  b.87 

5.37  b.93 

b.O 

5.b7  b.oq 

5.66  7.00 

5.0 

3.72  6.08 

12.9b  b8.03 

6.0 

11.38  23. bU 

109.78  232.50 

7.0 

bo.  13  62,58| 

bOl.07  625.66 

The  number  of  dose  levels  is  2k+l. 
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Table  6. 2. a 

Information  (I  )  for  Estimation  of  the  Mean  of  the  Normal  Tolerance 

'  a 

Distribution  with  Scale  Parameter  Known,  for  the  Finite  Experiment 


Middle 

Dose 

Location 

Dose 

Range 

N 

<r2I  (x  ) 
a  o' 

X  -Li 
o  ^ 

<*  ; 

k-1-^ 

k-2 

k-3 

k*U 

k»8 

k«20 

k*h0 

0 

1/2 

1.67 

1.6U 

1.63 

1.63 

1.62 

1.62 

1.62 

0 

1 

1.98 

1.87 

1.83 

1.82 

1.79 

1.78 

1.77 

0 

2 

3.1U 

2.82 

2.66 

2.59 

2.U8 

2.U2 

2.38 

0 

h 

U.70 

9.56 

5.16 

U.98 

U.71 

h.Sh 

U.hU 

0 

10 

U.71 

7.85 

10.81 

12.22 

11.76 

11.35 

11.08 

The  number  of  dose  levels  is  2k+l. 
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i.e.  the  greater  the  number  of  dose  levels,  the  greater  the 
information . 

N 

Table  8.2.b  presents  the  values  of  /a^j  (x  )  for 

av  o' 

values  of  a,  k  and  N,  for  values  of  xq/ jj,. 

8 .  U  Efficiency  of  the  Spearman  Estimator  for  the  Finite  Experiment, 
Scale  Parameter  Known 

The  efficiency  — ^  of  tiio  Spearman  estimator  for  the  finite 
experiment  will  be  defined  as  the  r  atio  of  the  inverse  of  the 
information  to  the  mean  square  error  of  the  estimator: 

Vi,  (x  ) 

e  Cx>  — 2— 2 —  (8-W 

MSE  (x  I  x  ) 
a  I  o' 

Table  8. 3. a  presents  computational  results  for  the  normal 

distribution  for  x  *u  for  several  values  of  a  and  k. 

o 

Table  8.3 .b  presents  computational  results  for  the  normal 
distribution  for  several  values  of  a,  k  and  ^  for  values  of 

V 

8.f?  Efficiency  of  the  Spearman  Estimator  for  the  Finite  Experi¬ 
ment,  Scale  Parameter  Unknown 

Table  8.3.b  indicates  that  the  efficiency  of  the  Spearman 


1 

The  efficiency  (E  )  as  defined  in  terms  of  information 
(8.U)  is  a  useful  measure  because  the  information  is  intrinsic 
to  the  experiment  itself  and  not  dependent  on  any  method  of 
estimation.  The  inverse  of  the  information  cannot  be  taken 
as  an  absolute  lower  bound  on  the  variances  of  all  estimators, 
nor  can  it  be  assumed  that  there  is  any  estimator  with  vari¬ 
ance  ti.is  small.  However,  it  is  a  lower  bound  for  the  vari¬ 
ances  v f  all  unbiased  estimators. 
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Table  8.2.b 

Information  (I  )  for  Estimation  of  the  Mean  of  the  Normal  Tolerance 

cl 

Distribution  with  Scale  Parameter  Known,  for  the  Finite  Experiment 


Middle  1 

Dose 

Location 

Dose 

Range 

N 

a2I  (x  ) 
a  0' 

X  -u. 

(  °  ) 
cr  ' 

(±  -) 
o 

k*l-/ 

k*2 

k«a 

0 

±1 

1.98 

1.8? 

.5 

2.10 

2.00 

1.0 

2.U9 

2.U3 

2.0 

5.13 

5-53 

3-0 

20.51 

25.08 

iuO 

197.37 

270.27 

0 

±2 

3-3U 

2.82 

.5 

3-33 

2.86 

1.0 

3-36 

3.01 

2.0 

3.90 

U.09 

3.0 

6.62 

8.55 

h.C 

22.78 

314.18 

5.0 

; 

205. U8 

328.95 

0 

—b 

5.56 

a. 98 

.5 

5.5U 

a. 98 

1.0 

5.52 

a. 98 

2.0 

5.56 

5.02 

3.0 

5.61 

5.a2 

a.o 

1 

6.51 

7.37 

5.o 

11.03 

15.39 

6.0 

i 

37.97 

61.52 

7.0 

3U2.U7 

592.11 

The  number  of  dose  levels  is  2k +1. 
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estimator  can  be  very  small  for  M»100  when  the  a  priori 
estimate,  x  ,  for  p  is  in  error.  The  information  was  com- 
puted  assuming  scale  parameter  known.  If  the  scale  para¬ 
meter  is  unknown  the  information  is  considerably  decreased 
for  xQ/p.  Table  8.ii  presents  the  element  I^(x  )  of  the 

inverse  of  the  information  matrix  for  estimation  of  the 
mean  of  trie  normal  tolerance  distribution  for  both  location 
and  scale  parameter  unknown.  Table  8.$  presents  the  cor¬ 
responding  efficiency  (E^)  of  the  Spearman  estimator  for 

cl 

two  unknown  parameters.  Tables  6. 3. a,  8.3.b  and  8.5  demon¬ 
strate  that  for  the  uvaal  finite  level  design  with  limited 
numbers  of  subjects  the  Spearman  estimator  has  high  efficiency 
relative  to  the  information  in  the  experiMnt,  when  the  toler¬ 
ance  distribution  is  normal. 

9.  REGULAR  BEST  ASYMPTOTICALLY  NORMAL  ESTIMATORS  WITH  THE  WRONG 
MODEL 

9.1  General  Discussion 

One  advantage  of  the  Spearman  estimator  is  that  no  para¬ 
metric  form  need  be  specified  for  the  tolerance  distribution. 
This  advantage  would  be  of  no  practical  importance  if  the 
conpeting  parametric  estimator  based  on  a  common  model  has  a 
distribution  that  is  insensitive  to  moderate  changes  in  the 
functional  form  of  the  true  tolerance  distribution.  In  this 
section  computations  are  presented  to  illustrate  the  effect 
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Tabic  8. 3. a 

Efficiency  of  the  Spearman  Estimator  for  the  Mean  of  the  Normal 
Tolerance  Distribution  with  Scale  Parameter  known,  for  the  Finite 
Experiment 


Middle 

Dose 

Location 

Dose 

Range 

Ea(x  )-/ 
a  0 

X  -u  J 

<-§-> 

(± 

k-1^/ 

k-2 

k-3 

k-U 

k*8 

k-20 

k-ljO 

0 

1/2 

3.27 

14.56 

5.09 

5. U3 

6.00 

6.U8 

6.75 

0 

1 

1.28 

1.58 

1.73 

1.82 

1.97 

2  .OU 

2.13 

0 

2 

.9h 

1.00 

1.03 

1.03 

1.05 

1.08 

1.09 

0 

h 

.39 

.9h 

1.02 

.99 

.99 

.93 

.99 

0 

10 

.06 

.25 

.56 

.80 

.95 

.99 

.99 

y  E  (x  ) - i - 1 - 

I  (x  )MSE  (x  x  ) 
a  o  a  *  o' 

2/ 


The  number  of  dose  levels  is  2k+l. 
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Table  8.3.b 

Efficiency  (E  )  of  the  Spearman  Estimator  for  the  Mean  of  the  Normal 
Tolerance  Distribution  with  Scale  Parameter  Known,  for  the  Finite 


Experiment 


Middle 

Dose 

Location 

"X  -p 

(-S-) 

! 

Dose 

Range 

(±  -) 
a 

E  (x 

a' 

o)y 

k-i-/ 

N-10 

k-2 

k-1 

! 

N-100 

k-2 

k*U 

o 

±1 

1.28 

1.58 

1.28 

1.58 

•5  i 

1.39 

1.61* 

1.09 

.85 

1.0 

1.61 

1.1*9 

.55 

.28 

2.0 

1.00 

.69 

.11 

.07 

3.0 

.88 

.80 

.09 

.08 

U .  o 

3.15 

3.57 

.32 

.36 

0 

±2  i 

•  9k 

1.00 

1 

.91* 

1.00 

o 

; 

i 

.99 

1.02 

.99 

1.02 

1.0 

! 

l.Ol* 

1.11 

1 

l.ol* 

1.09 

2.0 

1.19 

1.73 

1.12 

.76 

3.0 

! 

2.51 

1.57 

.56 

.18 

U.o 

1 

2.03 

1.U6 

! 

1 

.21 

.15 

5.0 

I 

5.12 

5.26 

.51 

.53 

0 

±1* 

.91* 

.98 

.,9k 

.98 

.5 

.98 

.98  j 

.98 

.98 

1.0 

1.02 

,98 

) 

1.02 

.98 

2.0 

.9k 

.99 

.9k 

.99 

3.0 

l.Ol* 

1.11 

l.ol* 

1.10 

U.o 

i 

1.19 

1.81* 

1.15 

1.05 

5.0 

i 

2.97 

2.53 

.85 

.32 

6.0 

3.3U 

2.62 

.36 

.26 

7.0 

i 

8.53 

9.U6 

.85 

.95 

y 


W 


i _ 

\>0 


2/ 

The  number  of  dose  levels  is  2k  +1* 


Table  8.1* 


Inverse  Information  (i^)  for  Estimation  of  the  Mean  of  the  Normal 

A 

Tolerance  Distribution  with  Scale  Parameter  Unknown,  for  the  Finite 


Experiment 


Kiddle 

Dose 

Location 

) 

Dose 

Range 

!  (±  |) 

i 

ran(x  )y 

a  o 

k-1-^ 

k-2 

k-a 

0 

±1 

1.98 

1.87 

.5 

2.1*0 

2.53 

1.0 

j 

U. 31 

5.a3 

2.0 

! 

38.63 

51.20 

3-0 

1 

907.36 

959.53 

0 

±2 

3.3U 

2.82 

.5 

3.33 

2.87 

1.0 

3.37 

3-13 

2.0 

1 

a .  7i 

6.98 

3.0  j 

66.80 

62.9a 

0 

±U 

5.56 

a. 98 

.5 

5.5a 

a. 98 

1.0 

5.52 

a.  98 

2.0 

5.56 

5.03 

3.0 

5.62 

5.63 

U.o 

7.85 

12.57 

5.0 

III.33 

113.30 

„  (x  )  is  the  element  of  the  inverse  of  the  information 
a  '  o ' 


y  iu(x 

a  < 

matrix  corresponding  to  tne  estimator  of  [i  . 
The  number  of  dose  levels  is  (2k+l). 
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Table  8.5 

Efficiency  (E*1)  of  the  Spearman  Estimator  for  the  Mean  of  the  Normal 
& 

Tolerance  Distribution  with  Scale  Parameter  Unknown,  for  the  Finite 
Experiment 


Middle 

Dose 

Location 

X  -M- 

<:  j 

Dose 

Range 

E11 

a 

y 

<V 

N-10 

k*l“/  k»2  k-U 

N-100 

k-1  k»2  k»l* 

0 

±1 

1.28  1.58 

1.28  1.58 

1.59  2.08 

1.25  .88 

1.0 

2.79  3.3U 

.96  .39 

2.0 

7.52  6. ill 

.82  .66 

3.0 

38 .9k  30.76 

3.90  3-08 

0 

±2 

.91*  1.00 

.91*  1.00 

.5 

.99  1.03 

.99  1.03 

1.0 

1.05  1.15 

1.05  1.12 

2.0 

1.1*3  2.95 

1.35  2.53 

3.0 

25.29  26.56 

5.63  11.55 

0 

±1*  . 

.91*  .98 

.91*  .98 

.5 

.98  .98 

.98  .98 

1.0 

1.02  .98 

1.02  .98 

2.0 

.91*  .99 

.91*  .99 

3.0 

1.05  1.16 

1.05  1.1 1* 

i*.o 

1.1*1*  3.H* 

1.39  1.79 

5.o 

29.93  18.65 

8.61  2.36 

y  -n 


E~~(x  ) 
a  N  o/ 


,11,  X 

\  (xo} 
MSEa(x  |  xQ) 


If 


The  number  of  dose  levels  is  2k  +1 
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on  some  RBAN  estimators^-  due  to  changes  in  the  functional 
form  of  the  tolerance  distribution. 

The  angular  model  is  used  and  the  characteristics  of  some 
RBAN  estimators  based  on  this  model  are  examined  for  true  tol¬ 
erance  distributions  with  the  forms:  logistic,  normal  and 
uniform.  The  experimental  designs  for  which  computations  are 
given  are  one,  two,  and  five  level  designs  with  scale  para¬ 
meter  known,  and  some  two  level  designs  with  scale  parameter 
unknown. 

The  RBAN  estimator  used  in  each  of  these  finite  designs 
is  an  explicit  function  of  the  independent  binomial  variates 
corresponding  to  the  several  dose  levels.  Consequently  the 
mean  and  variance  of  the  limiting  normal  distribution  of  the 
estimator  (as  the  sample  sizes  at  the  fixed  dose  levels  in¬ 
crease)  can  be  computed.  These  values  are  called  the  asymp¬ 
totic  mean  and  asymptotic  variant"'.  Since  the  estimator  is 
inconsistent  when  the  wrong  model  is  used,  the  asymptotic 
mean  square  error  is  computed  from  the  asyirptotic  mean  and 
variance  and  this  value  is  compared  with  the  asymptotic  vari¬ 
ance  of  the  RBAN  estimator  under  the  correct  model. 

9 m2  One  Level  Experiment 

The  model  used  for  the  tolerance  distribution  is  the 
angular  distribution  (see  Table  2.2).  Assume  that  the  scale 


See  Ncyman  (26)  for  the  definition  of  RBAN  (regular 
best  asyrptotically  normal)  estimators*  See  Taylor  (29)  for 
a  discussion  of  RBAN  estimators  in  bio&ssay. 
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parameter,  3  *  **  known  and  that  the  experiment  consists  of 

v 

testing  N  subjects  at  the  dose  level,  x-0.  Let  the  observed 

proportion  responding  at  x-0  be  denoted  by  p  and  the  expected 

—1  ,  — 

proportion  be  denoted  by  P.  Let  y*sin~  Jp  -  M/^  .  Then  the 
maximum  likelihood  estimator  of  |j.  is: 


|i*  -  (9.1) 

1  P 

The  asymptotic  mean,  variance  and  mean  square  error  of  the 
estimator,  for  a  given  value  of  P,  are: 


VaA(4>  ’^2  (CN) 
HSEaA(P>  -p 


where  Y-sin 


-1 


(9.2) 

(9.3) 

(9.U) 


and  P-E(p) . 

If  the  true  tolerance  distribution  is  angular  with  scale 
parameter  ^  ,  then  p-  -  ~  and  the  asymptotic  mean  square 
error  is  the  variance  (9.3).  Denote  the  angular  distribution 
by  G.  If  the  true  tolerance  distribution  is  F/G,  then  the 
estimator  will  not  be  consistent.  The  asymptotic  variance 
will  remain  the  same  but  the  bias  contribution  to  the  mean 
square  error  will  not  be  zero. 
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The  asymptotic  mean  square  error  can  be  computed  for  any 
given  F  and  N.  Denote  this  value  by  MSE^Cp  j  F).  In 
specifying  F  it  is  necessary  to  choose  the  value  of  the  scale 
parameter.  This  should  be  done  so  that  the  F  is  "comparable" 
to  the  model,  G,  with  its  known,  fixed  ^  .  Three  methods 
for  choosing  scale  parameters  for  the  tolerance  distributions 
are  used:  (i)  equating  standard  deviations,  (ii)  equating 
the  distances  between  two  specified  percentiles,  (iii)  equat¬ 
ing  the  information  per  observation. 

As  a  measure  of  the  effect  of  the  tolerance  distribution, 
F,  on  the  estimator  based  on  G,  the  asymptotic  efficiency, 

E  .(G|F),  is  computed 

MSE  (n*  I  G) 

F)  -  - — -  . 

KSEaA(H*l  F) 

Tables  9.1. a,  9.1.b  and  9.1.c  contain  computational 
results  for  F  normal,  logistic  and  uniform  when  the  model  is 
angular.  The  results  indicate  that  the  differences  due  to 
the  several  tolerance  distributions  are  negligible  even  for 
N  of  100,  when  the  tolerance  distributions  are  equated  on 
distance  between  the  20°  and  6C  '  percentiles  and  the  ex¬ 
pected  proportion  responding  is  not  too  far  from  $0  percent, 
say  between  1 5  percent  and  ^5  percent.  However,  when  the 
distributions  are  equated  on  standard  deviation  the  asymptotic 
mean  square  error  does  show  a  marked  decrease  due  to  bias  for 
N  of  100  for  some  values  of  P  (e.g.  for  the  logistic,  E*.66 
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Table  9.1. a 

Asymptotic  Efficiency  (E  . )  for  an  RBAN  Angular  Estimator  far  Three 

&A 

Tolerance  Distributions  having  Standard  Deviations  the  Same  as  the 
Model  (One  Dose  Level  with  N  Subjects) 


Expected  Per 
Cent  Respond¬ 
ing  at  the  : 
Single  Dose  1 

EaA(G|F)i/ 

N  -  10 

N  «  100 

Level 

Logistic 

(F) 

Normal 

(F) 

Uniform 

(F) 

Logistic 

(F) 

Normal 

(F) 

Uniform 

(F) 

* 

1 

.1*3 

.66 

.69 

.07 

.17 

.18 

3 

.92 

.96 

.89 

.56 

.71 

.1*5 

5 

1.00 

1.00 

.97 

.99 

1.00 

.78 

7 

.96 

.99 

1.00 

.7U 

.93 

.97 

10 

.91 

.97 

1.00 

.51 

.79 

.96 

15 

.87 

.96 

.97 

.1*0 

.69 

.78 

20 

.87 

.95 

.96 

.3U 

.68 

.69 

25 

.89 

.95 

.96 

.1*5 

.72 

.68 

30 

.92 

.97 

.96 

.51* 

.78 

.72 

35 

.95 

.98 

.97 

.66 

.85 

.79 

ao 

.98 

.99 

.99 

.81 

.93 

.88 

1*5 

.99 

1.00 

1.00 

.91* 

.98 

.97 

5o 

1.00 

1.00 

1.00 

1 

1,00 

1.00 

1.00 

1/ 


W°lF) 


i 

1  F) 


where  F  Is  the  true  tolerance  distribution  andG  is  the 
angular  model.  See  (9.1)  for  the  definition  of  **^  . 
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Table  9.1.b 

Asymptotic  Efficiency  (E^)  for  an  RBAN  Angular  Estimator  for  Three 
Tolerance  Distributions  having  Distances  Between  the  20th  and  60th 
Percentiles  the  Same  as  the  Model  (One  Dose  Level  with  N  Subjects) 


Expected  Per  I 
Cent  Respond*] 


ing  at  the 
Single  Dose 

N  -  10 

| 

r* - - 

N  -  100 

Level 

tiogistic 

(F) 

Normal 

(F) 

Uniform  j 

(F) 

Logistic 

(F) 

Normal 

(F) 

Uniform 

(F) 

1 

.15 

.37 

.50 

.02 

.06 

.09 

3 

.Uo 

.63 

.68 

.06 

.18 

.18 

5 

.62 

.8U 

.81 

.1U 

.3U 

.30 

7 

.78 

.92 

.88 

.27 

.53 

.UU 

10 

.92 

.97 

.95 

.53 

.78 

.68 

15 

.99 

1.00 

.99 

.92 

.97 

•9h 

20 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

\r\ 

CNJ 

.99 

1.00 

11.00 

.98 

.99 

.99 

30 

.99 

1.00 

1.00 

.96 

.99 

.97 

35 

1.00 

1.00 

1.00 

.97 

1.00 

.97 

Uo 

1.00 

1.00 

1.00 

.99 

1.00 

.99 

U5 

1.00 

1.00 

1.00 

.99 

1.00 

1.00 

50 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1/ 

\a(0  1  F) 

I'SE  , 

« 

I'SE 

SJi 

o£|o> 

where  F  is  the  true  tolerance  distribution  and  G  is 

# 

the  angular  model.  See  (9.1)  for  the  definition  of  n^. 


Table  9»I*c 


Asymptotic  Efficiency  (E^)  for  an  RBAN  Angular  Estimator  for  Three 
Tolerance  Distributions  having  Information  the  Sane  as  the  Model 
(One  Dose  Level  with  N  Subjects) 


Expected  Per 
Cent  Respond¬ 
ing  at  the 
Single  Dose 
level 

Z 

U(QlF)i/ 

II  -  10 

1 

N  -  100 

Logistic 

(F) 

Normal 

(F) 

Uniform 

(F) 

Logistic  Normal 
(F)  (F) 

Uhiform 

(F) 

1 

.11 

.01 

.01 

.01 

.00 

.00 

3 

.20 

.oh 

.01* 

.02 

.00 

.00 

5 

.30 

.09 

.10 

.oh 

.01 

.01 

7 

.50 

.17 

.19  | 

.07 

.02 

.02 

10 

.58 

.32 

.36  j 

.12 

.Oh 

.06 

15 

.61* 

.6*1 

.70 

.15 

.15 

.19 

20 

.80 

.71 

.75 

.28 

.20 

.23 

25 

.89 

.86 

.88 

• 

4^ 

VA 

.38 

.1*2 

30 

.95 

.9k 

.95 

.66 

.&) 

.61* 

35 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

hO 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

h5 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

50 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1/ 

E  .(GIF 

vni? 

)  -  ^ 

(4f°> 

£La  ' 

f  aA 

COl  t  F) 

where  F  is  the  true  tolerance  distribution  and  G  is 
the  angular  model.  See  (9.1)  for  the  definition  of 
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9.3 


at  P«3$  percent,  and  E—.3U  at  P»20  percent).  Equating  in¬ 
formation  shows  even  larger  effects  on  the  efficiency  for 
snail  values  of  P. 

The  results  for  the  angular  model  and  estimator  have 
been  corroborated  by  repeating  the  confutations  for  the  log¬ 
istic  model  and  its  estimator.  Results  analogous  to  those 
in  Table  9.1.b  are  presented  in  Table  9.2  for  the  logistic 
estimator. 

Since  the  results  presented  are  asymptotic  approxima¬ 
tions,  it  is  of  interest  to  see  whether  the  relationships 
indicated  by  these  computations  are  valid  in  the  range  of 
sample  size  used.  Exact  confutations  analogous  to  Table  9.1. c 
are  presented  in  Table  9.3  for  N=10. 

Two  Level  Experiment 

Again  let  the  model  and  estimator  be  based  on  the  angular 
distribution  (Table  2.2).  Let  the  experiment  consist  of  f!J 
subjects  tested  at  x*x^*  and  subjects  tested  at 

.  Let  p..  and  be  the  observed  and  expected  proportions 


respectively,  i*i,2.  Let  yu »sin”x  J"p^  - 
If  the  scale  parameter,  ^  ,  is  known, 


h  * 

the  RBAN  estimator 


used  will  be  denoted  by  and 


is: 


^2 " 

where  y  - 


Table  9.2 


Asymptotic  Efficiency  (E^)  of  an  RBAN  Logistic  Estimator  for  Three 
Tolerance  Distributions  having  Distances  Between  the  20th  and  the 
80th  Percentiles  the  Sane  as  the  hodel  (One  Dose  Level  with  N 
Subjects ) 


Expected  Per 
Cent  Respond¬ 
ing  at  the 
Single  Dose 
Level 

W°  1  F> 

17 

r - 

N  -  10 

N  -  100 

formal 

(F) 

Angular 

(F) 

Uhi/ara 

(F) 

Normal 

(F) 

Angular 

(F) 

thlfozm 

(F) 

1 

.91) 

.79 

.65 

.63 

.27 

.16 

3 

.96 

.83 

.67 

.70 

.32 

.17 

5 

.97 

.88 

.71) 

.79 

Ji2 

.22 

7 

•99 

.92 

.81 

.86 

.55 

.30 

10 

.99 

.96 

.90 

.93 

.71) 

.1)8 

15 

1.00 

1.00 

.90 

.99 

.9U 

.85 

20 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

25 

1.00 

1.00 

1.00 

1.00 

.98 

.91) 

30 

1.00 

1.00 

1.00 

1.00 

.97 

.89 

35 

1.00 

1,00 

1.00 

1.00 

.97 

38 

1)0 

1.00 

1.00 

1.00 

1.00 

.98 

.93 

1)5 

1.00 

1.00 

1.00 

1.00 

.99 

.98 

50 

1.00 

1.00 

1.00 

1.00 

1.00 

X00 

3^(0  |  F) 


I  °> 

|  F) 


where  F  is  the  true  tolerance  distribution  and  G  la 
the  logistic  nodel. 


75 


Table  9.3 


Exact  Efficiency  (E^)  for  an  REAM  Angular  Estimator  for  Three  Tol¬ 
erance  Distributions  having  Information  Equal  to  that  for  the  Model 
(One  Dose  level  with  10  Subjects) 


Expected  Per 
Cent  Respond¬ 
ing  at  the 
Single  Dose 
Level 


1 

3 

5 

7 

UO 

15 

20 

£5 

30 

35 

Uo 

U5 

50 


.06 

•17 

.29 

•U1 

.56 

.76 

.89 

.95 

.98 

.99 

1.00 

1.00 

1.00 

]ISEa(n£  |  G) 
MSEfi(^  |  F) 


.01 

.06 

.18 

.37 

.66 

•9U 

1.01 

1.01 

1.01 

1.00 

1.00 

1.00 

1.00 


.01 

.06 

.21 

.uu 

.76 

.98 

1.02 

1.02 

1.01 

1.00 

1.00 

1.00 

1.00 


where  F  is  the  true  uclerance  distribution  and  G  is 
the  angular  nodel. 


If  the  scale  parameter  is  not  known,  the  RBAN  estimator 


used  will  be  denoted  fcy  ^  and  is: 


y2"yl 


(9.6) 


The  asymptotic  means,  variances  and  mean  square  errors 
for  the  two  estimators,  for  given  and  P ^  are: 


Y  +Y 

EaA(4}  "  ‘  — 

*  2  2(i 


(9.7) 


where  I 


i-  ^_1/^  A 


(9.8) 


Y  +Y 

X1  2x2 


(9.9) 


EaA^2  ^ 


VY2 

2(Y2-Yi) 


(9.10) 


2  2 
Y  +Y 
2  *1 

(v  -y  ')  ^ 
u2  l' 


(9.11) 


2  2 
Y  +7 

*' x  1  X2  X1 


USE  (li,.,  )  •  *-T 
ajr^2  ;  2h 


<w 


♦ 


Y  +Y 
2  *1 

2(Y2-Yx) 


(9.12) 


If  the  t  olorance  distribution  is  angular,  ^  is  consist¬ 


ent!  and  if,  in  addition,  the  scale  parameter 


is  corsistent. 


is  p 


then  \x.r 
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Denote  the  angular  tolerance  distribution  by  0.  If  the 
true  tolerance  distribution  is  FjA},  then  both  estimators  will 
be  inconsistent.  The  asymptotic  mean  square  error  can  be 
determined  for  given  F  and  N,  for  specified  values  of  and 

£ 

P^.  Tablos  $lln-a  and  9.U.b  present  asymptotic  efficiencies  of  ^ 

■ic  t 

and  for  F  taken  to  be  logistic  when  the  model  is  angular. 

The  results  in  Tables  9.U.a  and  9.U.b  show  that  in  both  cases 

the  effect  of  a  change  in  F  from  angular  to  logistic  has 

little  effect  on  the  asymptotic  mean  square  error  of  the 

angular  estimator. 

9-h  The  2k +1  Level  Experiment 

Again  let  the  model  and  estimator  be  based  on  the  angular 

N 

distribution.  Let  the  experiment  consist  of  subjects 

tested  at  each  of  2k+l  levels.  Let  the  dose  levels  be 

1  -k+1  1  0  1  k-1  1 

"  r  :k  ’  *  *  *  -  2k-'  2k’  •  *  *  Ik  ■ •  2  * 

Let  the  observed  proportions  be  denoted  by  p^  and  the  expected 

proportions  by  P^.  Let  y^«  sin  1  Jp^  -  . 

If  the  scale  parameter,  ^  ,  is  known,  then  the  maximum 

* 

likelihood  estimator  is  denoted  by  and  is 


*  y 

Ur'  *  - 

5  p 


where  y 


k 

-k  3 

2&T 


(?. 13) 


The  asyr^ptotic  mean,  Variance  arid  mean  square  error  of  the 
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Table  9.h.a 

Asymptotic  Efficiency  (E^)  for  an  RBAW  Angular  Estimator  for  the 


Logistic  Tolerance  Distribution  (Scale  ,  Parameter  Known-:  Two  Dose 
Levels,  Subjects  at  each  Dose  Level) 


Expected  Per  Cent  Responding 

(G  \t)l/ 

Level 

<V 

Level 

(P2) 

N  -  10 

N  »  100 

10 

t 

1$ 

.70 

.19 

10 

20 

*80 

.28 

10 

30 

.93 

.56 

10 

ho 

.98 

.88 

10 

5o 

1.00 

1.00 

10 

60 

.99 

*9h 

10 

70 

.98 

.88 

10 

80 

.98 

.89 

10 

90 

1.00 

1.00 

20 

ho 

1.00 

.97 

20 

60 

1.00 

.99 

20 

80 

1.00 

j 

1.00 

ho 

50 

1.00 

1.00 

ho 

60 

1.00 

1.00 

where  F  is  the  logistic  tolerance  distribution  and  G 


is  the  angular  model.  See  (9-5)  for  the  definition 
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Table  9.k.b 

Asymptotic  Efficiency  (E^)  for  an  RBAN  Angular  Estimator  for  the 
Logistic  Tolerance  Distribution  (  Scale  Parameter  Uhknown:  Two 
Dose  Levels,  Subjects  at  each  Dose  Level  ) 


Expected  Per  Cent  Responding; 
Lower  Upper  j 

Level  Level 

(Pl>  <P2> 

EaA(G 

N  -  10 

N  ■  100 

10 

; 

1.00 

.96 

10 

20 

•99 

.92 

10 

30 

.99 

.91 

15 

ko 

1.00 

.96 

10 

50 

1.00 

1.00 

IP 

60 

.99 

.96 

10 

70 

.98 

.89 

10 

80 

.98 

.89 

10 

90 

1.00 

1.00 

20 

kO 

1.00 

.99 

cO 

60 

1.00 

1.00 

20 

80 

1.00 

1.00 

20 

kO 

50 

1.00 

1.00 

ko 

60 

1.00 

1,00 

y 


2^(0  S  F) 


1  G) 

^4'  i « 


where  F  is  the  logistic  tolerance  distribution 
and  G  is  the  angular  model*  Tee  (9.6)  for  the 

•t 

definition  of  . 
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estimator  are,  for  given  P^: 

V#  •  *  T  ^ 

where  Y^»sin*’^  J~P^  - 
k 

-k  Yi 

and  Y  ■  - 

2k +1 


VaA^> 


(9.15) 


MSE  (\£)  -  — ♦  (p*  ~)2  (9.16) 

**  *  UNp*  P 

When  the  distribution  is  angular,  the  estimator  is 
consistent.  Let  the  angular  tolerance  distribution  be  denoted 
by  G.  For  F/G,  the  asymptotic  mean  square  error  can  be  calcu¬ 
lated.  For  G  located  to  give  a  specified  set  of  a  com¬ 
parable  F  must  be  chosen.  An  F  cannot  be  chosen  which  will 
give  the  same  values  as  the  model  at  all  levels.  In  the  com¬ 
putations  F  was  chosen  to  give  the  same  values  at  the  end¬ 
points  of  the  range  of  dose  levels,  i.e.  at  and  x^. 

Table  9.5  presents  the  results  of  confutations  for  the  case 
of  five  dose  levels,  angular  estimator,  and  the  logistic  tol¬ 
erance  distribution.  The  effect  of  the  F  on  the  asymptotic 
mean  square  error  is  again  seen  to  be  negligible.  In  fact  the 
results  for  five  levels,  ^  known,  duplicate  almost  exactly  the 
results  for  two  levels, ft  known  (Table  9.1*.a). 
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Table  9.5 

Asymptotic  Efficiency  (E^)  ef  an  RBAK  Angular  Estimator  for  the 
Logistic  Tolerance  Distribution  (  Scale  Parameter  Known:  Five  Dose 
Levels,  N/5  Subjects  at  each  Dose  Level  ) 


Expected  Per  Cent  Responding 
Lowest  Highest 

Dose  Dose 

<V 

E^(0  |F)3/ 

N  -  10 

N  •  100 

10 

' 

15 

i 

.70 

.19 

10 

20 

.79 

.27 

10 

30 

.91 

.50 

10 

hP 

.96 

.7  U 

10 

50 

.99 

.92 

10 

60 

1.00 

.99 

10 

70 

1.00 

1.00 

10 

80 

1.00 

1.00 

10 

90 

1.00 

1.00 

y  E^to !  r) 


tCEoA(M-S  I  0) 
*e,a(^;  i  D 


where  F  is  the  logistic  tolerance  distribution 

and  G  is  the  angular  model.  See  (9*1*0  for  the 

* 

definition  of  (i,  . 
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9.5  Summary  of  Wrong  Model  Investigation 

The  above  computations  indicate  that  in  the  case  of  para¬ 
metric  estimation  a  model  can  be  chosen  and  slight  deviations 
of  the  true  tolerance  distribution  from  the  nssinjiil  functional 
form  will  have  little  effect  on  the  mean  square  error  of  the 
estimator.  (It  is  important  that  the  tolerance  distribution 
have  about  the  same  spread,  measured  in  interpercentile  de¬ 
viation,  anticipated  in  the  model,  if  the  scale  parameter  is 
assumed  known.) 

Previous  sections  (7>  8)  have  shown  that  the  Spearman 
estimator  is  a  very  efficient  estimator  conpared  with  the 
amount  of  information  available  in  the  experiment.  This 
section  indicates  that  this  efficiency  would  not  be  greatly 
increased  if  comparisons  were  based  on  the  mean  s  quare  error 
of  fully  efficient  parametric  estimators  taking  into  consider¬ 
ation  the  possibility  of  using  a  wrong  model.  In  spite  of  this 
apparent  robustness  of  the  parametric  estimators,  the 
Spearman  estimator  la  recaasended  for  use  in  most  quantal  assay 
experiments  because  of  its  simplicity  and  high  efficiency. 
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APPENDIX  I.  THE  DISTRIBUTION  FUNCTION  AID  INFORMATION  FORIHE 
INFINITE  EXPERIMENT 

1.1  Distribution  Function  for  the  Infinite  Experiment 

Leitana  1.1  Let  F(x)  be  a  tolerance  distribution.  Specify  a  value 
for  d.  Choose  xq  randomly  from  the  interval  (0,d).  Let 

x  =x  +id,  i»0.±l,±2, .  Take  n  observations  on  a  Bernoulli 

variable  with  expected  value  F(x^)  for  i»0,±l, ±2, . . . .  Let 
all  the  observations  be  independent.  Let  r^  be  the  number 
of  responses  at  x^.  Then  this  infinite  experiment  determines 
a  probability  function  for  the  on  the  infinite  sample 
space. 

Proof:  Since  the  r^  are  independent  sums  of  independent 
Bernoulli  variables,  the  distribution  functions  for  finite 
sets  of  t  he  r^  satisfy  the  consistency  conditions  on  page  2 9 
of  Kolmogorov  (2h).  The  lemma  follows  irrsnediately  from  the 
theorem  on  the  same  page. 

1.2  Information  for  the  Infinite  Experiment,  Scale  Parameter  Known 

The  information  for  the  infinite  experiment  will  be 
defined  as  the  limit  as  k  — >  •  of  the  information  for  the 
finite  experiment,  with  dose  levels  xi«xQ+id,  i-O,  ±1,  ±2, 

...  ±k,  and  xQ  randomly  chosen. 

The  information  for  the  finite  experiment  with  f  ixed 

Xo  ia 

W  ■ E 

where  h  is  the  density  of  the  finite  sequence  of  binomial 
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variables  i-0,  ±1,  ±2,  ...  ±k. 


r. 


h  -  TT  (?i>  f/  d-Fi)n-ri 

-k  ^ 


k  nF 


I,  (x  )  -  £  vhere  T  .  *  -X  I 

k  0  -k  r1(i-Fi)  *  7v  x**i 

The  information  for  the  infinite  experiment,  I,  will 
be  defined  as 


I  •  lim  E 
k — >•  *o 


I.  (x  ) 
k'  o' 

0 


1  •  lim 
k  —*« 


,  k  nF  . 

i  2  ■ 

d 


dx 


-k  Fid-Pi) 


F 

- 1—  dt 

F(l-F) 


It  can  be  shown  that  if  n  is  a  translation 

parameter  then  I  is  finite. 

3.3  Information  for  the  Infinite  Experiment,  Scale  Parameter 
Unknown 

Let  F  be  written  F  ^  (x-n)J  .  The  information  matrix 
for  the  infinite  experiment  will  be  d  efined  as  the  matrix 
obtained  as  k— »•  for  the  finite  experiment. 
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where 


A 


1 


,'2 


F(l-F) 


dt 


F12 


J  F(l-F) 


dt 


(t-A1)2 


,'2 


F(l-F) 


nit 


f  F’2 

— -  dt 

J  F(l-F) 
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AFPENDIX  II.  BERNOULLI  PERIODIC  FUNCTIONS  AND  THE  2ULER- 
MACLAURLf  FORMULAE 

The  Bernoulli  Periodic  Functions 

The  Bernoulli  periodic  functions,  Pn(t),  and  the  Euler- 
HacLaurin  formulae  are  presented  here.  For  details  of  the 
development  of  the  f  ormulae  the  reader  is  referred  t  o  Cra¬ 
mer  (11)  pp.  122-125. 

The  Bernoulli  periodic  functions  are: 


(t) 


2  cos  t  _ 

k«l  22n"1(kTT)2n 


n*l, 2)090 


P?n+l(t)*  2  - - -  n-0,  2>*** 

w  22n(kn)2ntl 

These  functions  satisfy 

p'  (t)  »  (-1)31"1  P  ,(t)  . 
n  n-1 

The  first  three  Bernoulli  functions  are  (see  reference  1, 

p,  138): 

t  0<t<l 

p2(t)  '  +  T  0<t<1 

t  t2  t3 

p^(t)  *  -  $■*  ♦g-  °<  t<  1 

In  paragraph  U.2.7  the  sup  P  (t)  is  desired  for  n»l,2  and  3- 

t 


Using  the  polynomial  expressions,  the  following  are  obtained 
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sup  P,(t)  -  \ 
t 

sup  P2(t)  - 

sup  ?-(t)  ■  .0080... 
t  J 


II. 2  The  Euler -Mac  La  ur  in  Formulae 

Let  x  and  d  be  constants  and  let  the  1  term  in  a 
o 

finite  sum  be  g(xQ+id).  Let  g  be  continuous  with  a  continu- 

t 

ous  derivative  g  .  Then  the  first  Euler -KacLaurin  formula  is: 


k 

2  g(x  +id ) 
-k  0 


k 

-j  g(xQ+xd)dx  +  |g(xQ-kd)  +  ^(x^kd) 
-k 


k 

-d  |  P1(x)g' (xQ+xd)dX 

-k 


If  g  has  continuous  derivatives  of  higher  orders  repeated 
integration  by  parts  gives  the  Euler-KacLaurin  formulae: 

g(xo-id)dx  +  |  g(xQ-kd)  +  |  g(xQ+kd) 


2  g(x  +id) 
-k 


-k 


s  B 

-2 
1 


21 


,2i-l 


(2i)  / 


g(2l"1)(xo-kd)-g(2i"1)(xo>kd) 

P2s+1U)g(2s’1)(x0+xd)dx 


where  s  is  any  non-negative  integer,  must  exist  for 
j*l,2, . . „ (2s+l),  and  B..  are  the  Bernoulli  numbers,  defined  by: 
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The  first  Euler-MacLaurin  formula  glvea  above  can  be  applied 
to  an  infinite  sum  to  obtain: 

f  f 

Z  g(xo+id)  *  \  g(xQ+xd)dx  -  d  \  P^(x)g  (xQ+xd)dx 

_a»  »« 

provided  that  the  series  and  the  two  integrals  converge. 
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APPENDIX  III.  ESTIMATION  OF  THE  VARIANCE  CF  THE  SPEARMAN  ESTIMATOR 
III.l  An  estimator  for  V(x) 

An  obvious  estimator  for  the  variance  of  the  Spearman 
estimator  is: 

s2  »  2  p.q.  (if  converge*) « 

£  n-1  —  1  1  £ 

Then 


The  f  ollowing  results  can  b  e  obtained  for  this  estimator 
for  the  infinite  experiment  with  random  choice  of  x  ,  follow¬ 
ing  the  same  methods  used  for  obtaining  the  characteristics 
of  x, 

2 

If  F  has  a  first  moment,  s-  converges  with  probability 
one  and  has  the  following  mean  and  mean  square  error. 


E(s  f) 

x 


F(l-F)dx 


V(x) 


MSE(s2) 

x 


(1)3 

n 


2 

n  -6n+6 
n(n-l) 


F2(l-F)2dx 


F(l-F)dx 

+  Ex 

’b2(4|*0>’ 

J  j 

0 

2  2  2 
where  B‘(s^  J  xq)  denotes  the  bias  of  s^  as  an 

estimate  of  7(x)  (5*U),  conditional  on  x  . 

Unlike  the  case  of  x,  where  it  can  be  shown  that  the 

2 

MSE(x)  is  a  minimum  for  n«l  when  n'»n/d  is  fixed,  for  s- 


9U 


the  optimum  choice  of  n  for  n'  fixed  will  depend  on  F.  For 

2 

the  normal  tolerance  distribution  MSE(s  decreases  as  n 
increases  but  the  decrease  is  negligible  for  n  greater  than 
1*.  In  this  case,  then,  the  optimum  designs  for  estimation 
of  p  and  of  tf(x)  do  not  agree,  but  &  good  compromise  would 
be  to  choose  n-ii,  say. 

III. 2  An  Alternative  Estimator  for  ?(x)  Based  on  the  Second 
Moment  of  F 

Table  5*1  shows  V(x)  to  be  a  function  of  the  scale 
parameter : 

V(x)  -  ~  Cp  (III.1) 

where  Cp  is  a  constant  depending  on  the  parametric  formula¬ 
tion  of  the  tolerance  distribution.  The  constants  for  the 
normal,  logistic,  angular  and  uniform  distributions  are 
.561*2,  .5513,  .5750  and  .577U  respectively. 

If  the  constant  Cp  is  considered  over  all  distribution 
functions,  the  function  defined  ty : 

1  2 

P(x  ■  •■*■■■*  a  >  0 

a  1+a2 

P(x  *  a)  * 

1+a 

will  have  CL,-  — ,  so  that  C_  is  arbitrarily  close  to  zero 
F  lnT  F 

There  is  a  unique  maximising  function  for  Cp 


for  small  a. 
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over  any  given  finite  interval  (see  Rustaggi/,  and  it  can 

be  shown  by  a  calculus  of  variations  argument  that  the 

uniform  distribution  is  this  function  (Cp».577ii). 

*2  _ 

An  alternative  estimator  (s-  )  for  7(x)  suggested  by 
(III.l)  is: 


(Vd/2-x)2 


^pi+l"pi^ 


The  estimate  of  a  used  in  s—  is  the  Spearman  type  of  esti¬ 
mator  suggested  'ey  Epstein  and  Churchman  (Hi)  and  shown  by 
Cornfield  and  Mantel  to  be  an  algebraic  approximation  to  the 
maximum  likelihood  estimator  of  cr  for  F  logistic. 


See  Rustagl#  Jagdish  Sharan,  "On  Minimizing  and  Maxi¬ 
mizing  a  Certain  Integral  with  Statistical  Implications," 
Annals  of  Mathematical  Statistics,  28  (1957)  309-328. 
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